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SUMMARY 

Algorithms  and  procedures  are  presented  for  the  reduction  of  force  and  moment  data 
derived  from  wind-tunnel  models  supported  by  internal  strain-gauge  balances.  The 
algorithms  are  developed  in  their  most  general  forms,  suitable  for  implementation  on  the 
new  generation  of  powerful  minicomputers  currently  being  included  in  wind-tunnel  data 
acquisition  systems.  Although  the  emphasis  of  the  report  is  on  the  treatment  of  data  derived 
from  sting  mounted  strain-gauge  force  balances,  the  analysis  is  applicable,  with  only 
minor  modifications,  to  data  derived  from  modern  load-cell  based  external  force  balances. 
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arbitrary  vector 

coefficient  in  transformation  matrix  ABB. 
speed  of  sound 
slot  width 

coefficient  of  transformation  matrix  ABB. 

tunnel  width 

span 

cross  wind  force 

tunnel  cross-sectional  area 

sensitivity  and  interaction  coefficient 

sensitivity  and  interaction  coefficient  matrices 

chord 

drag  force 

diagonal  matrix  of  primary  sensitivity  coefficients 

control  surface  length  scale 

slot  spacing 

slot  parameter 

arbitrary  force  vector 

fuselage  finesc  ratio 

arbitrary  couple 

unit  gravity  vector 

total  pressure 

tunnel  height 

total  applied  load 

component  load 
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h 

hinge  moment 

[J] 

critical  stress  coefficient  matrix 

!*■] 

deflection  matrix 

ki 

element  of  [K] 

L 

lift  force 

1 

rolling  moment 

1 

length 

M 

Mach  number 

m 

pitching  moment 

n 

jawing  moment 

P 

static  pressure 

P 

porosity  parameter 

<? 

dynamic  pressure  (=  \/2pv‘l) 

R 

strain  gauge  amplifier  output 

Re 

Reynolds  number 

S 

reference  area 

T 

interference  parameter 

l 

aerofoil  thickness 

u 

velocity  component 

V 

volume 

V 

free  stream  velocity 

V 

velocity  component 

W 

weight 

M' 

velocity  component 

X 

axial  force 

[n  [ni,  [X2] 

normalized  interaction  coefficient  matrices 
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ordinate 

Y 

side  force 
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normal  force 
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P 
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4> 
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angle  of  incidence 
angle  of  sideslip 

Prandtl-Glauert  compressibility  parameter  =  (I  —M2)* 

ratio  of  specific  heats 

interference  parameter 

approximation  to  applied  load  H 

blockage  factor 

accelerometer  outputs 

angular  deflection 

pitch  angle 

direction  cosine  matrix 
dynamic  viscosity 
kinematic  viscosity 

angular  rotation  from  balance  to  body  axes 
density 

balance  critical  stress 
interference  parameter 
angular  deflection 
roll  angle 
angular  deflection 
yaw  angle 

angular  rotation  from  accelerometer  to  body  axes 


Superscripts 


non-linear  vector  or  matrix  of  component  loads 
transpose  of  a  matrix 


(  F 

[  1' 


matrix  inverse 


Subscripts 


accelerometer  axes 
body  axes 
balance  axes 
base 

body  or  fuselage 
free  air 
hinge 

missile  axes 

wind  off  value 

rig  or  indicated  value 

stability  axes 

sine  definition 

solid  body 

tailplane 

tunnel  axes 

tangent  definition 

uncorrected  (i.e.  measured)  values 

wind  axes 

wake 


Note — Symbols  designating  forces  and  moments — X,  Y,  Z.  I,  m,  «,  D,  C,  L  and  h  are  also  used 
as  subscripts  when  forces  and  moments  are  expressed  in  coefficient  form.  e.g.  C i.  -  lift 
coefficient. 
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I.  INTRODUCTION 


An  operating  wind  tunnel  is  potentially  capable  of  generating  large  quantities  of  raw  data 
in  a  very  short  period.  Hie  ever  increasing  cost  of  staffing,  maintaining  and  running  major  wind 
tunnels  has  generated  a  strong  economic  inducement  to  increase  wind  tunnel  efficiency  b> 
speeding  up  the  acquisition  and  recording  of  test  data.  One  common  method  of  achieving  this  aim 
has  heen  to  use  computer  based  data  acquisition  and  processing  systems.  From  the  time  when 
general  purpose  digital  computers  first  became  available,  they  have  been  used  in  such  systems, 
l.arly  (late  I950's  and  early  1960's  I  computer  based  systems  made  use  of  off-line  general  purpose 
machines.1  These  machines  were  generally  remote  from  the  wind  tunnel,  data  being  recorded 
at  the  wind  tunnel  site  (usually  with  relatively  slow  paper  tape  or  card  punches)  and  transported 
to  the  central  machine  for  later  processing.  With  such  systems,  results  were  often  not  available 
for  hours  or  even  days  after  the  test. 

The  advent  of  the  minicomputer  in  the  mid-1960  s  led  to  the  development  of  dedicated 
wind  tunnel  computer  systems.  Processed  lest  results  were  then  available  immediately  following 
the  test,  and  in  some  cases  in  real-time  (or  nearly  so)  during  the  lesi.  Connection  of  the  computer 
directly  to  the  wind  tunnel  measuring  instruments  (often  using  parallel  data  transfers)  also  allowed 
a  significant  increase  in  the  rale  at  which  data  could  be  recorded  and  processed.  Large  increases 
in  tunnel  efficiency  resulted  from  the  increased  data  collection  rate,  hut  more  importantly 
from  the  new-found  ability  to  detect  and  correct  erroneous  data  while  a  test  was  in  progress. 
However,  carfy  dedicated  minicomputer  systems  suffered  from  slow  execution  speeds  and  from 
severely  limited  memory  capacity,  (his  led  to  the  need  to  minimize  the  size  of  data  reduction 
programs  and  sometimes  to  undesirable  limitations  on  generality  or  compromises  on  overall 
accuracy. 

Developments  in  the  computer  technology  since  the  1960  s  have  seen  the  cost  effectiveness 
of  minicomputers  increase  by  factors  of  up  to  50.  In  current  installations  it  is  not  uncommon 
I  r  a  dedicated  wind  tunnel  minicomputer  to  support  4  Megabytes  of  memory  and  to  operate 
at  speeds  of  one  million  operations  per  second  or  more.  This  class  of  machine,  which  is  signi¬ 
ficantly  more  powerful  than  even  the  largest  of  the  machines  available  in  the  1960's,  allows  a 
complexly  general  approach  to  be  taken  when  implementing  data  acquisition  and  processing 
algorithms.  No  longer  there  any  excuse  for  computer  generated  compromises  on  accuracy. 
These  new  machines  also  allow  system  designers  to  explore  quite  new  approaches  to  wind  tunnel 
data  acquisition.  For  example,  the  high  speeds  of  these  machines  permit  a  much  greater  degree 
of  real  time  data  monitoring  and  active  model  attitude  control,  leading  to  the  feasibility  of 
completely  new  methods  of  acquiring  data,  such  as  constant  parameter  testing  techniques. - 
Indeed,  it  is  now  feasible  for  a  wind  tunnel  test  to  be  conducted  completely  under  computer 
control 

The  purpose  of  this  report  is  to  set  out.  in  the  most  general  way.  the  algorithms  and  pro¬ 
cedures  for  the  acquisition  and  processing  of  wind  tunnel  data  derived  from  internal  strain 
gauge  force  balances.  Although  some  material  has  been  directly  reproduced  from  other  sources, 
it  was  considered  worthwhile  to  bring  together,  in  their  most  general  forms,  all  the  algorithms 
required.  With  the  more  powerful  minicomputers  now  in  use  in  wind  tunnels,  such  an  approach 
may  be  implemented  almost  directly.  The  development  assumes  that  the  major  source  of  data 
is  the  familiar  internal,  sting  mounted  strain  gauge  balance,  but  the  analysis  could  be  applied, 
with  only  minor  modifications,  to  data  acquired  from  modern  load-cell  based  external  balances. 

The  report  is  split  into  two  parts.  The  first  part,  consisting  of  Sections  2  and  3.  presents 
formal  definitions  of  the  axes  systems  in  common  use.  and  the  sets  of  angles  which  relate  one 
axes  system  with  another.  The  use  of  these  angles  to  transform  data  from  one  axes  system  to 
another  are  also  included.  The  second  part,  making  up  the  remainder  of  the  report,  sets  out  the 
procedures  and  algorithms  required  to  convert  raw  wind  tunnel  data  to  a  form  suitable  for  use 


I 


by  a  customer.  The  three  sections  which  make  up  this  part  of  the  report,  diside  the  overall 
process  into  three  logically  separate  stages : 

(i)  Calculation  of  the  aerodynamic  loads  acting  on  the  model.  This  involves  the  conversion 
of  strain  gauge  bridge  outputs  to  applied  loads,  correction  of  these  applied  loads  for 
the  effects  of  model  and  balance  mass,  transformation  to  body  axes,  and  formation  of 
non-dimensional  coefficients.  Allowance  must  also  be  made  for  the  effects  of  base 
pressures. 

(ii)  Determination  of  the  attitude  of  the  model  and  hence  the  orientation  of  the  body  axes 
system  with  respect  to  the  relative  vsind  sect  >r, 

(iii)  Application  of  collections  for  the  effects  of  wind  tunnel  wall  interference  (carried  out 
in  the  wind  axes  system)  and  finally  the  transformation  of  data  to  the  axes  systems 
required  by  the  customer. 

Although  the  overall  approach  presented  in  this  report  is  followed  by  the  computer  based 
data  acquisition  system  recently  installed  in  the  A.R.L.  low-speed  wind  tunnel.3  this  system 
does  not  as  yet  fully  implement  all  of  the  procedures.  It  is  hoped  that  this  full  implementation 
can  be  achieved  in  the  near  futuic. 


2.  SYSTEMS  OE  AXES 

When  formulating  and  solving  problems  involving  internal  strain  gauge  balances  and  their 
mounting  systems,  a  number  of  frames  of  reference  (systems  of  axes)  must  be  considered.  Each 
of  the  axes  systems  considered  here  is  an  orthogonal  right-handed  system.  In  all.  six  a-es  systems 
will  be  considered.  This  is  not  an  exhaustive  set.  the  choice  being  based  on  convenience,  useful¬ 
ness.  and  to  some  extent  personal  preference. 

Tunnel  Axis.  vt.  it-  Jt  :  with  the  vt  axis  coincident  with,  and  in  the  same  direction  as  the 
relative  free-stream  wind  vector  (usually  horizontal),  and  the  it  axis  perpendicular  to  the  longi¬ 
tudinal  tunnel  axis  and  horizontal,  positive  to  the  right  when  looking  in  the  direction  of  .vt. 

Tunnel  axes  are  fixed  in  the  wind  tunnel  and  serve  as  the  locally  applicable  earth-fixed  (and 
inertial)  reference  frame.  Tunnel  axes  are  rarely  used  to  express  aerodynamic  loads. 

Wind  Axes.  t  w.  .  u.  :  with  the  vw  axis  parallel  to  the  relative  wind  vector,  positive 
into  the  wind,  and  the  :#  axis  in  the  plane  of  symmetry  of  the  model  if  it  has  one.  otherwise  arbi¬ 
trary. 

Wind  axes  are  the  traditional  frame  of  reference  for  the  expression  of  aerodynamic  coefficients 
and  the  calculation  of  wind  tunnel  wall  interference  (although  tunnel  axes  may  be  preferable  for 
the  latter  iT  the  tunnel  is  of  non-circular  cross-section  and  the  model  plane  of  symmetry  is  not 
parallel  to  one  of  the  tunnel  axes  of  symmetry).  Note  that  tunnel  axes  and  wind  axes  differ  by 
at  most  an  angle  of  roll.  If  the  model  does  not  have  a  plane  of  symmetry,  it  is  customary  to  take 
the  Y w  axis  to  be  horizontal  when  the  model  is  in  its  defined  zero  roll  orientation. 

Body  Axes.-— xh.  in.  m:  with  the  vu  axis  along  the  model  longitudinal  axis,  positive  for¬ 
ward.  and  the  zn  axis  in  the  model  plane  of  symmetry,  positive  "downwards''. 

Body  axes  remain  fixed  in  the  model  and  remain  related  to  it  via  reference  planes,  e.g.  the 
chord  or  fuselage  reference  line,  plane  of  symmetry  and  plane  of  the  wings  for  an  aircraft  model. 
For  models  without  planes  of  symmetry,  e.g.  missiles,  the  i  n  and  zn  axes  are  usually  defined  to 
be  in  some  convenient  mutually  perpendicular  directions. 
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Stability  Axes. — xs,  >%  zs:  with  the  xs  axis  parallel  to  the  projection  of  the  relative  wind 
vector  on  the  plane  containing  the  xb  and  >>b  axes,  positive  forward,  and  the  >>s  axis  coincident 
with  and  in  the  same  direction  as  the  yB  axis. 

Stability  axes  are  a  special  set  of  body  axes  of  considerable  importance  in  the  study  of  flight 
dynamics.  The  stability  axes  yaw  with  the  model  but  do  not  pitch. 

Missile  Axes. — Xm.  >m,  zm:  with  the  xm  axis  coincidant  with  and  in  the  same  direction  as 
the  ,vb  axis,  and  the  ym  axis  coincident  with  and  in  the  same  direction  as  the  jb  axis  with  the  yaw 
and  roll  angles  zero. 

Missile  axes,  also  known  as  "non-rolling  body  axes",  are  frequently  used  in  missile  flight 
dynamics. 

The  origins  for  all  of  the  above  axes  systems  are  usually  assumed  to  coincide  when  the  model 
is  in  its  reference  (zero  pitch,  yaw  and  roll)  condition.  The  position  of  the  origin  in  this  reference 
condition  is  normally  chosen  to  be  the  location  of  the  full-scale  aircraft  centre  of  gravity  in  the 
model,  but  may  be  any  convenient  reference  position.  Depending  on  the  details  of  the  equip¬ 
ment  used  to  change  model  attitude,  the  origins  of  the  various  axes  systems  may  no  longer 
coincide  at  finite  model  attitudes. 

In  general,  the  moment  centre  used  for  the  calibration  of  a  strain  gauge  balance  will  not 
coincide  with  the  chosen  origin  for  the  body  axes  system.  In  order  to  allow  use  of  balance  cali¬ 
brations  carried  out  in  arbitrary  axes  systems,  it  is  useful  to  define  a  further  set  of  axes  - 
balance  axes. 

Balance  Axes.--x' b,  .1  n,  z'b  :  w  ith  the  x’b  axis  parallel  to  the  longitudinal  axis  of  the  balance 
or  other  axis  used  for  the  x-direction  during  calibration,  positive  "forwards",  and  the  z'b  axis 
similarly  defined,  positive  "downwards". 

The  balance  axes  system  is  equivalent  to  the  body  axes  system,  but  with  the  origin  at  the 
moment  centre  assumed  during  balance  calibration  and  with  the  axes  in  those  directions  used 
during  balance  calibration.  Note  that  it  is  usual  to  assume  that  the  balance  is  infinitely  stiff  between 
the  “model  end"  and  the  moment  centre,  i.e.  the  origin  and  direction  of  the  balance  axes  are 
fixed  with  respect  to  the  model  end  of  the  balance. 

For  all  the  above  axes  systems,  the  components  of  the  net  force  F  acting  through  the  origin 
are  defined  to  be  X.  Y  and  Z  in  the  x,  y  and  z  directions  respectively  (with  suitable  subscripts 
added  to  indicate  the  particular  axes  system  in  use).  For  compatibility  with  aircraft  stability 
theory,  three  additional  forces  L  (lift).  D  (drag)  and  C  (cross-wind  force)  are  defined  in  the  wind 
axes  system  (and  occasionally  in  the  stability  axes  system).  These  additional  forces  are  related  to 
X.  Y  and  Z  by 


L  -  Z  w 
D  =  Aw 
C  =  -  Fw 

The  components  of  the  net  couple  G  acting  about  the  origin  are  defined  as  I.  m  and  n. 

Note  that  in  the  remainder  of  this  report,  the  subscript  B  (body  axes)  will  be  dropped  w  here 
there  is  no  possibility  of  confusion. 

In  the  body  axes  system,  the  components  of  the  relative  wind  vector  V  in  the  vh.  An  and  zb 
directions  are  defined  to  be  u,  r  and  >v  respectively. 


3.  DEFINITIONS  OF  ANGLES 

Reduction  of  wind-tunnel  data  involves  the  use  of  several  frequently  used  angles.  These 
angles  fall  into  two  categories  :  aerodynamic  angles  and  positioning  or  orientation  angles. 
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3.1  Aerodynamic  Angles 

The  orientation  of  the  model  with  respect  to  the  relative  wind  vector  V  can  be  defined  in 
terms  of  either  the  three  orthogonal  components  of  V  in  the  body-axis  system  ( u ,  r.  «  ),  or  more 
conveniently,  by  the  magnitude  of  V  and  two  suitably  defined  angles.  These  two  angles  are  known 
as  the  angles  of  incidence  x,  and  sideslip  fi.  Currently  there  are  two  alternate  definitions  in  use 
for  x  and  fi  -the  tangent  definitions  and  the  sine  definitions — and  there  is  a  need  for  precision 
in  stating  which  definition  is  in  force. 


Tangent  Definitions 

xi  is  the  angle  between  the  model  .VB-axis  and  the  projection  of  the  wind  vector  on  the  xb. 
-  h  plane,  i.e. 

iv 

tan  xt  — 

a 

fit  is  the  angle  between  the  model  ,vB-axis  and  the  projection  of  the  wind  vector  on  the  xb. 
>b  plane,  i.e. 

tan  fit  =  V 

ti 


Sine  Definitions 

x»  is  the  angle  between  the  wind  vector  and  its  projection  on  the  model  xb.  .vb  plane,  i.e. 
tv 

s.n  «s  =  y 

fis  is  the  angle  between  the  wind  vector  and  its  projection  on  the  model  xb.  cr  plane,  i.e. 


sin  /8S  = 


r 

V 


The  merits  of  the  two  sets  of  angles  thus  defined  (and  combinations  thereof)  have  been 
considered  by  Warren.1  As  he  points  out,  the  traditional  approach  to  aircraft  stability  and  control 
was  to  use  the  wind  axes  system,  and  to  assume  that  squares  of  component  velocities  in  the  t  B 
and  zn  directions  (r2,  iv2)  could  be  neglected,  and  that  the  v-component  of  velocity  u  was  equiva¬ 
lent  to  the  total  velocity  V.  Under  these  conditions,  tangent  and  sine  definitions  are  equivalent. 
Since  that  time,  testing  of  missiles  and  modern  aircraft  at  large  combined  values  of  x  and  fi 
(i.e.  large  values  of  v  and  »  )  has  invalidated  these  assumptions  and  hence  the  two  definitions  must 
be  differentiated. 

Workers  in  the  missile  flight  dynamics  field  prefer  to  use  the  same  definition  of  both  inci¬ 
dence  and  sideslip  to  preserve  the  inherent  symmetry  of  the  missile  configuration.  However 
the  choice  of  either  the  tangent  or  sine  definition  is  rather  arbitrary,  although  there  is  some  argu¬ 
ment  for  using  the  tangent  definition.  In  the  aircraft  field,  the  increasing  use  of  stability  axes 
in  stability  and  control  theory  leads  to  the  preference  for  the  tangent  definition  of  incidence  and 
the  sine  definition  of  sideslip.  This  choice  (xt,  (S„)  also  leads  to  particularly  simple  definitions  in 
terms  of  the  orientation  angles  for  the  traditional  yaw /pitch  rig  commonly  used  in  low-speed 
wind  tunnels  (see  Appendix  A). 
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For  the  remainder  of  this  report  relationships  will  be  developed  only  for  the  combination 
on,  ds  and  their  subscripts  will  be  dropped  where  there  is  no  possibility  for  confusion.  However, 
relationships  for  other  combinations  will  be  presented  in  the  Appendices  so  that  that  combination 
best  suited  to  a  particular  test  (or  customer)  may  be  used. 


3.2  Orientation  Angles 

The  orientation  of  one  set  of  axes  relative  to  any  other  set  can  be  given  by  three  angles,  which 
are  the  consecutive  rotations  about  the  axes  z,  >\  x  in  that  order  that  carry  one  set  of  axes  into 
coincidence  with  the  other.  This  is  a  special  case  of  “Euler  Angles”.  In  wind  tunnel  testing  the 
Euler  angles  of  interest  are  those  that  rotate  the  tunnel-axes  system  xt,  it.  -t  into  coincidence 
with  the  body-axis  system,  and  the  angles  are  </>,  0  and  4>.  Referring  to  Figure  I,  these  angles  may 
be  defined  as : 

0  -the  yaw  angle — rotating  the  tunnel-axes  system  about  ozt  to  give  the  axes  .va,  .ra,  za- 

W— the  pitch  angie— rotating  the  x>,  >’a,  za  system  about  oyi  to  give  the  axes  X3.  .13,  Z3. 

4>  -the  roll  angle — rotating  the  X3,  >3,  Z3  system  about  0.V3  to  give  the  body-axes  system 
a  h.  yn,  zb- 

The  positive  directions  for  t/i,  dand  <t>  obey  the  normal  right-hand  screw  rule  about  the  relevant 

axes. 

It  must  be  noted  that  the  sequence  in  which  the  rotations  (>/>.  0,  <j>)  are  taken  is  important. 
If  the  sequence  is  changed,  the  orientation  of  the  body-axes  system  with  respect  to  the  (fixed) 
tunnel-axes  system  (and  hence  the  values  of  a  and  /?)  will  also  change. 

To  avoid  ambiguities  which  can  otherwise  result  in  the  set  of  angles  (|4,  (>.  $)  the  ranges  are 
limited  to 


or  0  ^  1(1  <  2  n 
it/2  ^  tt  it  2 


or  0  ^  <t>  $  2tt 


The  Euler  angles  arc  then  unique  for  most  model  orientations  and  gi\e  aerodynamic  angles 
in  the  ranges 


n  ^  j  £  n 


and 


tt  n 


3.3  Relationships  Between  Aerodynamic  and  Orientation  Angles 

In  the  wind  tunnel  the  model,  and  hence  the  body-axes  system  can  he  rotated  through  the 
Euler  angles  (</i,  0,  6),  and  this  rotation  causes  the  aerodynamic  angles  1  and  rt  to  exist  As  pointed 
out  by  Warren,4  perhaps  the  clearest  way  of  showing  the  angles  of  incidence  and  sideslip  that  are 
produced  by  the  rotation  (0,  0,  4>)  is  to  visualize  the  trace  of  a  point  on  the  tn-axis  on  the  surface 
of  a  sphere  whose  centre  is  at  the  origin. 
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Figure  2  shows  the  general  case.  If  P  is  the  initial  position  of  the  point  on  the  .va-axis,  the 
.vb,  )’b  plane  being  the  horizontal  plane  WPE,  and  the  .vb.  zb  plane  being  the  vertical  plane  NPS, 
rotating  the  model  through  an  angle  of  yaw  p  moves  P  to  Pi,  the  .vb,  >b  plane  remaining  in 
WPPiE  and  the  ,n»,  zb  plane  moving  to  NP\S.  Rotating  through  an  angle  of  pitch  0  moves  P\ 
to  Pi,  the  .vb,  >b  plane  moving  to  W\P,E\  while  the  .vb,  zb  plane  remains  in  NP>P\S.  Finally 
rotating  through  an  angle  of  roll  p  does  not  change  Pi,  but  the  .vb.  vb  plane  moves  to  WiPiEt 
and  the  .vb,  zb  plane  to  P-±N->S->.  If  PQ  and  PR  are  perpendicular  to  KiPiS-i  and  HiP-iEi  respec¬ 
tively,  then  the  definitions  of  the  angles  of  sideslip  and  incidence  are  equivalent  to  arcs  of  great 
circles  on  the  surface  of  the  sphere  as  follows: 

m  PiQ 


A  P-iR 


*s  =-  PR 

^  PQ 

By  the  application  of  elementary  spherical  trigonometry  it  may  be  shown  that 


sin  «s  -  cos  P  sin  0  cos  P- l-sin  p  sin  P 
sin  0*  =  cos  </«  sin  ^  sin  oi  sin  p  cos  p 
tan  p  sin  p 


tan  at  tan  0  cos  p  + 

tan  j3t  =  tan  fl  sin  P 


cos  0 

tan  p  cos  P 
cos  6 


(3.1) 

(3.2) 

(3.3) 


(3.4) 


Special  forms  of  the  above  equations  applicable  to  particular  cases,  together  with  their 
forms  in  the  presence  of  additional  small  deflections,  are  set  out  in  Appendix  A. 


3.4  Transformations  Between  Systems  of  Axes 


Consider  two  orthogonal  right-handed  systems  or  axes  denoted  V|.  x->,  V3  and  .Vi',  xi,  xp. 
The  direction  of  the  .ti',  v>',  tV  set  of  axes  with  respect  to  the  ,vi.  V2.  V3  set.  may  be  expressed  as  a 
matrix  of  direction  cosines. 


*11 

*12 

*13 

*21 

*22 

*23 

*31 

*32 

*33 

where  Ap(,  is  the  cosine  of  the  angle  between  the  v',,  and  v,|  axes. 

Amongst  the  many  useful  properties  of  this  matrix  which  may  be  derived,5  those  of  interest 
in  the  present  area  are  as  follows : 

(i)  The  inverse  of  the  direction  cosine  matrix  and  its  transpose  are  equal,  i.e. 

[A]  1  -  [A)T  (3.5) 
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(ii)  Consider  a  vector  A.  The  components  of  A  in  x\,  xz\  xz,  A'  where 


A.l 


A' 


A2’ 


A3' 


are  related  to  its  components  in  the  xi,  x-i,  xz  system.  A,  via 


A'  =  (A]A 


(3.6) 


This  transformation  may  be  readily  inverted  using  (3.1). 

A  =  [Aj-iA' 

-  MTA'  (3.7) 

In  the  present  context  we  are  interested  in  the  ability  to  express  the  model  loads  determined  in 
one  axes  system  in  some  other  axes  system.  For  example,  the  components  of  the  model  force 
vector  F  in  the  tunnel  axes  system  (Ft)  may  be  obtained  from  its  components  in  the  body  axes 
system  (FB)  from 


Fb  =  [A]bt  Ft 


(3.8) 


where  [A]Bt  is  the  matrix  of  direction  cosines  which  relate  the  direction  of  the  body  axes  system 
to  the  tunnel  axes  system.  It  may  be  shown  that  in  terms  of  the  Euler  angles  (1 p,  0,  </,)  which 
rotate  the  tunnel  axes  into  the  body  axes,  this  matrix  is  given  by 


'cos  0  COS  1/1 

cos  0  sin  >l> 

-sin  0 

sin  <P  sin  0  cos  1 P 

sin  4>  sin  0  sin  p 

sin  <P  cos  0 

—cos  <t>  sin  ip 

+  COS  <t>  COS  Ip 

cos  <f>  sin  6  cos  <P 

cos  <P  sin  0  sin  P 

cos  <P  cos  0 

_  +  sin  <£sin  i/j 

— sim£cos0 

The  inverse  transformation — from  tunnel  axes  to  body  axes  is  then  given  by 
Ft  =  [A]tb  Fb 
=  [A]bt"'  Fb 

=  [A]btt  Fb  (3.10) 

Appendix  B  presents  transformation  matrices  for  the  transformation  of  vector  components 
between  the  more  often  utilized  sets  of  axes — body,  stability  and  wind-axes.  The  matrices  are 
presented  for  all  useful  combinations  of  the  aerodynamic  angles  «»,  jS9,  «,  and  /St. 


4.  DETERMINATION  OF  AERODYNAMIC  LOADS  ON  THE  MODEL 

This  section  addresses  the  problem  of  the  determination  of  the  aero dynamic  loading  applied 
to  a  mode!  from  the  outputs  of  amplifiers  connected  to  each  strain  gauge  bridge  output.  It 
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will  be  assumed  in  what  follows  that  the  linearity  of  the  strain  gauge  bridge  amplifiers  is  sufficient 
for  them  not  to  effect  the  relationship  between  applied  load  and  amplifier  output,  and  hence  that 
a  balance  calibration  is  independent  of  the  particular  set  of  amplifiers  in  use  (except  possibly 
for  an  adjustment  for  amplifier  gain).  In  practice,  this  assumption  is  valid  for  the  amplifiers 
in  use  at  ARL,  the  newer  types6  being  an  order  of  magnitude  more  linear  (0  01  “0)  than  the 
expected  accuracy  of  a  balance  calibration  (0- 1  %). 

The  determination  of  aerodynamic  loads  may  be  considered  in  several  stages.  Firstly, 
experience  has  shown  that  for  the  newer  amplifiers,  although  the  overall  gain  of  the  bridge/ 
amplifier  system  is  extremely  stable,  the  system  is  susceptible  to  zero  drift  which  must  be  corrected. 
This  zero  drift  has  several  sources,  the  major  contribution  being  the  effect  of  temperature  on 
the  strain  gauge  bridge  itself  (both  bulk  temperature  and  temperature  gradients).  For  some 
older  amplifiers,  a  correction  for  amplifier  gain  drift  must  also  be  made. 

Secondly,  the  gross  loads  acting  on  the  model  must  be  calculated  by  the  solution  of  the  balance 
equations  for  a  six-component  balance  or  by  the  application  of  simple  linear  calibration  to 
hinge-moment  measuring  bridges.  Finally  the  model  weight  tares  must  be  calculated  and  sub¬ 
tracted  from  the  gross  loads  to  give  the  net  aerodynamic  loads  acting  on  the  model. 

The  components  of  aerodynamic  loads  thus  calculated  will  be  expressed  in  the  balance 
axes  system.  Before  being  useful  for  further  analysis,  the  load  components  must  be  converted 
to  dimensionless  coefficient  form,  and  the  influence  of  stings  on  pressures  acting  on  blunt  bases 
taken  into  account.  Finally  the  corrected  balance  axes  coefficients  must  be  transformed  to  the 
body  axes  system,  taking  into  account  both  the  angular  relationships  between  the  two  axes 
systems  and  any  offset  between  the  balance  moment  centre  and  the  origin  of  the  body  axes 
system. 


4.1  Corrections  to  Balance  Outputs 

For  a  model  mounted  on  a  six-component  strain  gauge  balance,  the  outputs  of  the  balance 
will  be  presented  as  six  output  channel  readings  Rt'  i  —  1.6.  If  hinge  moments  are  being  meas¬ 
ured  there  will  be  additional  outputs  Rh’  from  the  moment  measuring  strain-gauge  bridges-  - 
one  for  each  hinge  moment.  In  general,  these  outputs  will  not  be  zero  in  the  wind-olf  condition 
and  this  value  (the  “zero"  reading)  must  be  subtracted  from  the  wind-on  outputs 

Ri  ’  Ri  Rig  \ 

and  I  (4.1) 

Rh  Rh  Rito'  ) 


It  is  also  desirable  to  correct  the  raw  balance  outputs  for  any  drift  in  these  "zero"  readings 
during  a  run.  This  may  be  achieved  if  a  wind-off  "zero"  reading  is  recorded  at  the  same  altitude 
at  the  beginning  and  end  of  each  wind-tunnel  run  and  the  wind-on  outputs  corrected  via 


and 


Ri  -  Ri 


Rh  Rh’ 


[(Rig  IlH-RltmlnK+f Rio  lend  I 

2  J 

)beglnnlng  +  (  ^Ho  )onct  | 

2  J 


(4.2) 


Obviously  it  is  not  possible  to  apply  such  a  correction  scheme  until  the  run  has  been  com¬ 
pleted  and  the  wind-off  “zero"  reading  recorded.  For  data  presentation  in  “real-time"  (i.c. 
during  a  run),  the  correction  scheme  of  equation  (4.1)  must  be  used. 
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For  strain-gauge  readout  equipment  having  variable  sensitivity  and  where  a  "calibration" 
of  the  sensitivity  setting  is  possible,  a  scheme  analogous  to  that  of  (4.2)  can  be  applied  to  beginning 
and  end  of  run  "calibration  "readings. 


4.2  Reduction  of  Balance  Outputs  to  Gross  Forces  and  Moments 

Although  the  designer  of  a  strain-gauge  balance  attempts  to  make  each  output  sensitive 
to  only  one  load  component,  it  is  not  possible  to  eliminate  entirely  interactions  due  to  other 
components.  In  many  cases,  especially  where  the  balance  design  loads  are  large  for  the  size  of 
the  balance,  it  is  necessary  to  take  into  account  the  presence  of  non-linear  interactions  in  addition 
to  the  usual  linear  interactions. 

In  general  the  outputs  of  the  strain  gauge  bridge  for  each  component  are,  due  to  these 
interactions,  functions  of  all  component  loads.  The  reading  for  the  ith  balance  component  is 
usually  expressed  as  a  polynominal  function  of  the  applied  loads,  and  for  well  designed  balances 
it  is  generally  found  that  terms  of  third  and  higher  degree  in  load  components  are  negligible, 
while  the  coefficients  of  second  degree  terms  are  small,  but  possibly  significant.  Hence  each 
component  output  nv>y  be  represented  to  sufficient  accuracy  by  a  second  order  polynomial  of 
the  form 


Ri  —  Ci.iHi  +  C|  2W2  + . -f  Ci  sffe 

+  Cl,ll//l2  +  Cl,22H,22+ . +C\taHs2 

+  C'iii2//iF2  +  CijsHiHs  + . +Ci.MHsH$ 

Where  the  C's  are  the  calibration  coefficients  determined  during  balance  calibration  and 
H\  are  the  component  loads.  The  calibration  coefficients  may  be  classified  as  follows : 

(a)  "linear",  e.g.  Cij.  j  =  1,6 

(b)  "load  squared",  e.g.  Ci.u,  j  =  1,6 

(c)  “load  cross  product",  e.g.  Cij*,  j  =  1,5,  k  =  j+  1,6 

It  follows  that  for  each  component  there  are  27  calibration  coefficients,  giving  a  total  of  162  for 
a  six-component  balance. 

The  set  of  six  simultaneous  equations  represented  by  (4.3)  to  describe  balance  output  are 
most  conveniently  represented  in  matrix  form,  with  the  calibration  matrix  partitioned  into 
“linear"  and  “non-linear”  sub-matrices : 

R  =  [CIJ  H  +  [C2J  H*  (4  4) 


where 


R  is  the  vector  of  component  outputs 

fCl  j  is  the  matrix  of  linear  calibration  coefficients,  Cij 

[C2]  is  the  matrix  of  non-linear  calibration  coefficients,  Ci.jj  and 
Cl. Jk 

H  is  the  vector  of  component  loads 

H*  is  the  vector  of  squares  and  cross  products  of  component  loads. 
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T  * 


It  is  useful  to  "normalize"  the  linear  and  non-linear  calibration  matrices  [Cl  1  and  [C2]  by 
premultiplying  each  by  a  diagonal  matrix  (Oj  composed  of  the  diagonal  elements  of  the  linear 
calibration  matrix  [Cl  |  giving 


[All  [0|[CI| 
I A  2)  |/>J1C2] 


(4.5) 


This  has  the  effect  of  allowing  the  elements  of  [A  I)  and  |A'2)  to  have  universal  application  while 
those  of  |/>]  are  dependent  on  the  particular  readout  system  in  use.  A  factor  for  converting 
from  one  readout  system  to  another  of  the  same  type  can  then  be  simply  applied  to  the  “primary 
sensitivity"  matrix  (/>| 

The  procedure  lor  determining  the  component  load  vector  from  the  component  outputs  via 
equation  <4  4)  has  been  described  previously’  *-8  but  is  included  here  for  completeness.  Pre- 
multiplication  of  both  sides  of  (4  4)  by  [0|  yields 

[/>!»  [0||CI|H+(/)|[C2|H* 

or 

H  [AI|H»IA2|H* 

where  H  is  a  vector  of  approximate  component  loads  ignoring  all  (both  linear  and  non-linear) 
interactions.  Solving  for  H  gives 

H  =--  [A'l|  '  H  [AT  |  1  (A'2|  H*  (4.6) 

If  there  are  no  non-i,  .ar  terms  in  the  balance  calibration  then  [A  2|  0  and  H  is  given  by 

the  explicit  solution 


H  =  [A’l  1  '  H 


(4.7) 


If  however  non-linear  terms  are  significant,  then  (4  6)  must  be  solved  by  an  iterative  method. 
A  first  approximation  to  H  may  be  obtained  by  ignoring  the  non-linear  terms,  i.e.  from  (4.7) 

Hi  [Al|  1  H 

and  this  approximation  may  be  used  to  obtain  a  first  approximation  of  the  non-linear  term  (A,) 
thus 


A,  [VI  j  1  | A  2 1  H, * 

where  the  elements  of  the  vector  Hi*  are  composed  of  the  squares  and  cross-products  of  the 
elements  of  the  approximate  vector  Hi. 

A  second  approximation  Hz  is  then  obtained  from 

H2  -  H!  A, 

This  procedure  may  be  generalised  into  an  iterative  scheme  where  an  approximation  to  the 
non-linear  terms 


An  =  [AT]-'  [A’2J  Hn* 


(4.8A) 
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is  used  to  produce  a  new  approximation  to  the  component  loads  F  from 

Hn  —  Hi—  i)  (4.8B) 


where  n  is  iteration  number. 

The  convergence  of  this  scheme  has  been  investigated  in  Reference  7,  but  for  typical  balances 
with  small  non-linear  effects,  convergence  is  usually  very  rapid  (two  to  three  cycles). 

The  output  of  hinge  moment  measuring  bridges  can  usually  be  represented  to  sufficient 
accuracy  by  a  linear  calibration  of  the  form 


7?h  =  Ch  Hh 


(4.9) 


and  determination  of  the  applied  load  H h  from  the  output  R h  is  trivial. 


4.3  Tare  Weight  Corrections 

For  a  model  mounted  on  an  internal  strain-gauge  balance,  distinction  must  be  made  be¬ 
tween  the  part  of  the  applied  loads  determined  from  (4.8)  which  is  due  to  the  wind,  and  that  due 
to  a  change  in  model  attitude.  Due  to  the  finite  weight  of  a  model,  all  strain  gauge  balance  outputs 
will  include  the  combined  effect.  However,  since  the  attitude  of  the  model  with  respect  to  hori¬ 
zontal  is  known  or  can  be  determined,  the  balance  loads  due  to  the  weight  components  of  the 
model  can  be  calculated.  The  change  in  weight  components  due  to  the  change  in  attitude  from 
the  wind-off  “zero”  condition  to  the  wind-on  measurement  condition  may  then  be  calculated 
and  deducted  from  the  measured  (combined)  values  to  give  the  portion  due  to  the  wind  alone. 

The  distribution  of  the  model  weight  "tares”  among  the  balance  components  is  a  function 
of  the  magnitude  and  centre  of  gravity  position  of  the  “metric  mass"  and  the  balance  orientation 
with  respect  to  the  gravity  vector.  The  presence  of  terms  of  other  than  first  degree  in  the  balance 
calibration  matrix  introduces  a  further  difficulty  when  determining  tares.  This  is  the  problem  of 
defining  what  the  true  "zero”  of  the  load  scale  should  be.  For  anything  but  a  linear  calibration, 
the  use  of  one  "zero"  for  calibration  of  the  balance,  and  others  for  load  measurement  will  give 
rise  to  errors.  To  avoid  this  problem  the  concept  of  an  "absolute  zero  load"  or  “bouyant" 
condition,  first  introduced  by  Galway.9  will  be  used  when  calculating  the  tare  components. 
This  tare  condition  is,  by  definition  independent  of  the  balance  orientation  and  the  magnitude 
and  centre  of  gravity  position  of  the  metric  mass.  All  incremental  measurements  relative  to 
this  "bouyant"  tare  condition  will  then  be  proportional  to  load  on  an  absolute  scale. 

When  using  the  concept  of  “bouyant"  tares,  it  is  important  to  define  exactly  what  is  meant 
by  the  "metric  mass"  and  to  carefully  distinguish  between  the  "metric  mass"  and  the  weight 
of  the  model.  The  “metric  mass"  may  be  defined  as  ail  parts  of  the  modcl/balance  system  as  a 
whole,  whose  weight  results  in  outputs  on  any  balance  component,  the  magnitude  and  distri¬ 
bution  of  which  depend  on  the  orientation  of  the  model/ balance  system  relative  to  the  gravity 
vector.  It  is  probable  that,  due  to  the  construction  of  the  balance  and  the  placement  of  individual 
strain  gauge  bridges,  different  components  will  sense  different  "metric  mass"  contributions  from 
different  parts  of  the  balance,  and  hence  different  total  tare  weights  and  centre  of  gravity  locations. 

The  problem  of  determining  the  components  of  measured  outputs  which  are  due  to  the 
wind  is  therefore  reduced  to  determining  the  “bouyant"  condition  for  the  strain  gauge  balance. 
The  method  used  here  is  one  in  which  the  bouyant  condition  for  each  balance  component  is 
derived  from  the  wind-off  outputs  from  a  series  of  distinct  model/balance  orientations.  It  is 
therefore  unnecessary  to  determine  the  model  weight  and  centre  of  gravity  location  by  separate 
means. 

Consider  first  the  loads  seen  by  each  component  of  the  balance  due  to  the  weight  of  the 
metric  mass.  For  all  types  of  model  support,  the  only  angular  motions  of  the  model /balance 
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with  respect  to  the  gravity  vector  are  pitch  and  roll.  Hence  the  loads  due  to  the  weight  of  the 
metric  mass  will  be  given  by 


Fx  =  tf  x  sin  ft 
Fy  =  W y  cos  0  sin  4> 

Fz  —  Wz  cos  0  cos  <4 
G i  =  IVi(}i  cos  0  cos  6  — 2(  cos  S  sin  <4) 
Gm  =  —  iVmlZm  sin  0  +  .tm  cos  0  cos  <4) 
G„  =  tf'nlxn  cos  0  sin  <ji  +  I  n  sin  H) 


y  (4.10) 


where  the  (T’s  are  the  components  of  the  metric  mass,  and  the  .t,  y,  e's  are  the  locations  (measured 
in  the  balance  axes  system)  of  the  centres  of  gravity  of  these  components,  as  seen  by  each  balance 
component. 

For  a  pitch-yaw  type  of  model  support,  where  the  roll  angle  <j>  is  constrained  to  be  zero, 
these  loads  reduce  to 


Fx  =  -  H'x  sin  0 
Fy  =  0 

Fz  =  Wz  cos  ft 
Gi  =  Wi  yi  cos  0 
Gm  —  -  Wm  (Ztn  sin  B  +  Xm  cos  0) 
Gn  =  W„  >'n  sin  ft 


h  (4.11) 


Initially,  with  the  balance  and  model  set  (with  the  wind-off)  at  some  attitude  80,  <j>0  (usually, 
but  not  necessarily  80  =-  O',  <t>a  =  0C),  the  output  of  each  component  is  recorded  as  a  “zero" 
reading,  (/?i)0.  This  reading  is  not  the  output  the  balance  would  have  at  that  attitude  at 
“bouyant”  conditions,  but  it  may  be  considered  as  an  approximation  of  the  bouyant  condition 
with  all  tare  components  of  the  metric  mass  assumed  to  be  zero.  The  balance  and  model  are 
then  moved  to  a  different  attitude  (still  with  the  wind-off)  where  the  output  of  each  component 
is(fli),.  The  change  in  output 


Afi,  =  (/?,),  (/?,)„ 

may  then  be  used  to  calculate  the  change  in  applied  load  AH\  corresponding  to  the  change  in 
attitude  (using  the  calibration  matrix  and  the  iterative  scheme  (4.8)).  This  calculated  change  in 
applied  load  may  be  used  with  the  equations  of  (4. 10)  or  (4.11)  to  obtain  an  approximation 
to  the  components  of  the  metric  mass.  For  example,  from  (4.1 1)  the  change  in  the  x-component 
load  is  given  by 


AFx  =  Wy.  (sin  ft  sin  ft„)  (4. 1 2) 

and  hence  Wx  may  be  determined. 

These  approximate  components  of  the  metric  mass  may  in  turn  be  used  to  determine  a 
better  approximation  to  the  “bouyant"  outputs  of  each  balance  component.  Returning  to 
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the  ^-component  example,  having  computed  If'x  it  may  be  used  to  calculate  the  tare  load  com¬ 
ponent  at  fi  =  »0  using  (4.10)  thus 

Fx*  —  H'x  sin  t>„ 

and  once  all  the  Hi*  have  been  determined,  the  balance  calibration  matrix  can  be  used  to  give 
the  /?i*  using  (4.4).  The  change  in  output  A/t,  for  the  change  in  attitude  is  then  recalculated 
but  now  including  the  approximation  to  the  "bouyant"  tare  reading 


A«,  =  (/?,),  -Ki'  (4.13) 

which  may  be  used  to  calculate  a  more  accurate  estimate  of  the  components  of  the  metric  mass. 
The  whole  process  is  repeated  until  successive  approximations  to  the  metric  mass  components 
agree  to  within  some  tolerance. 

It  should  be  noted  that  the  above  procedure  does  not  yield  explicit  estimates  of  the  magni¬ 
tudes  of  all  the  metric  mass  components  or  of  their  centre  of  gravity  locations,  but  rather  esti¬ 
mates  of  the  values  of  their  products  (e.g.  fTmCm).  However,  since  these  products  are  in  fact 
the  values  required  in  the  application  of  the  overall  scheme,  this  does  not  represent  a  problem. 

In  principle,  only  two  model  balance  attitudes  are  necessary  to  determine  the  metric  mass 
components.  In  practice  however,  it  is  desirable  to  make  use  of  many  attitudes  to  increase  the 
accuracy  of  the  process.  The  approach  adopted  here  is  to  use  the  balance  outputs  from  N  distinct 
attitudes  and  then  to  represent  equations  such  as  (4.12)  as  least  squares  linear  approximations. 
Statistical  methods  may  then  be  used  to  calculate  standard  errors  of  estimate  for  each  component 
to  indicate  the  “goodness  of  fit”. 

Having  determined  the  required  values  of  the  components  of  the  metric  mass,  it  remains 
to  apply  these  and  the  concept  of  “bouyant"  tares  to  the  strain  gauge  balance  outputs  recorded 
during  a  wind  tunnel  run.  The  procedure  is  as  follows : 

(i)  The  loads  applied  to  the  balance  by  the  metric  mass  at  the  attitude  at  which  wind-off 
zeros  are  recorded  may  be  computed  from  equation  (4.10)  or  (4.1 1). 

(ii)  The  readings  of  the  six  balance  outputs  corresponding  to  these  loads  may  be  calculated 
from  the  balance  calibration  matrices  and  equation  (4.4). 

(iii)  Subtracting  these  readings  from  the  wind-off  zero  readings  /fi0'  gives  a  set  of  “bouyant" 
zero  readings. 

(iv)  These  bouyant  zeros  may  then  be  used  to  correct  wind-on  outputs  using  equation 
(4.1)  or  (4.2). 

(v)  The  resulting  set  of  absolute  readings  may  they  be  converted  to  a  set  of  equivalent 
absolute  loads  using  the  iterative  scheme  of  equation  (4.8). 

(vi)  Finally  the  aerodynamic  components  of  this  set  of  loads  is  calculated  by  subtracting 
the  tare  loads  due  to  the  metric  mass  at  each  wind-on  attitude  computed  from  either 
(4.10)  or  (4.11). 


it  should  be  noted  that  if  the  balance  calibration  matrix  contains  only  linear  terms,  i.e. 
[C 2]  in  equation  (4.4)  is  zero,  then  the  above  procedure  may  be  simplified.  In  this  case  it  is  suffi¬ 
cient  to  calculate  the  change  in  loads  on  the  balance  from  the  metric  mass  due  to  the  change  in 
attitude  from  the  wind-off  zero  to  the  wind-on  data  point.  This  change  in  load  is  then  deducted 
from  the  measured  wind-on  loads  to  give  the  portion  due  to  the  wind. 

The  output  of  hinge  moment  measuring  bridges  will  also  be  affected  by  the  “metric  mass" 
of  the  corresponding  control  surfaces.  For  most  conventional  control  surfaces  it  is  generally 
found  that  the  effects  of  the  metric  mass,  when  compared  to  the  required  accuracy  of  measure¬ 
ment  of  hinge  moments,  are  small  enough  to  be  ignored.  In  practice  the  wind-off  outputs  of 
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hinge  moment  measuring  bridges  should  be  observed  over  the  complete  range  of  attitudes  to 
be  tested.  If  the  magnitude  of  the  tare  effects  is  greater  than  allowed  by  the  desired  accuracy  of 
measurement,  an  analysis  similar  to  that  described  above  may  be  applied  to  determine  the  com¬ 
ponents  of  the  metric  mass  and  hence  to  separate  tare  and  wind  induced  calibrations  to  hinge 
moment.  Since  hinge  moment  calibrations  will  generally  be  linear  (e.g.  equation  (4.9)),  the 
simplifications  of  the  previous  paragraph  may  be  utilized.  Equations  similar  to  (4.10)  or  (4.1 1) 
giving  the  tare  contribution  to  hinge  moment  as  a  function  of  attitude  (in  this  case  including  the 
deflection  of  the  control  surface  itself)  may  be  formulated.  Special  attention  should  be  paid  to 
correctly  accounting  for  sign  conventions  for  both  hinge  moment  and  control  deflection  and  the 
effect  of  hinge  orientation  with  respect  to  the  overall  model  balance  axes  (e  g.  a  flap  or  aileron 
on  the  swept  trailing  edge  of  a  wing  or  a  model  mounted  on  a  pitch /roll  support  deserves  careful 
attention). 


4.4  Balance  Safe  Loads 

As  described  in  the  previous  section,  the  use  of  the  concept  of  "bouyant"  tares  allows  the 
calculation  of  the  absolute  loads,  i.e.  the  total  loads  resulting  from  all  external  fields,  acting  on 
the  balance.  Such  absolute  loads  are  exactly  those  required  for  checking  that  the  loading  state 
of  the  balance  is  such  that  maximum  stress  levels  remain  within  the  safe  operating  envelope. 
Traditionally,  it  has  been  the  practice  to  define  a  maximum  value  for  each  of  the  six  components 
applied  to  the  balance,  independent  of  the  magnitudes  of  the  other  five.  For  most  practical 
balances  this  will  be  a  very  conservative  approach  since  the  maximum  values  are  usually  derived 
on  the  basis  that  some  or  all  of  the  components  are  at  their  maximum  values  simultaneously. 
An  approach  which  allows  the  more  efficient  use  of  the  total  balance  capacity  is  to  check  that 
for  each  set  of  applied  loads  the  maximum  stresses  remain  within  safe  limits.  For  modern  bal¬ 
ances  for  which  no  on-balance  summing  of  bridge  outputs  is  used,  the  most  efficient  approach  is 
to  calculate  a  set  of  maximum  stresses  from  the  absolute  readings  produced  by  the  applied 
loads,  i.e. 


«  =  [/]R*  (4.14) 

where  o  is  a  vector  of  n  elements  corresponding  to  n  critical  stresses  and  \J]  is  an  nx6  matrix 
of  coefficients  determined  for  each  particular  balance  design,  a  is  then  compared  with  a  set  of 
n  maximum  allowable  stresses. 

This  approach  has  the  added  advantage  of  avoiding  the  time-consuming  iterative  solution 
required  to  determine  the  absolute  loads  H.  The  time  saved  is  specially  important  when  the 
state  of  the  balance  is  to  be  checked  at  regular  intervals  in  time  rather  than  only  at  the  times 
when  data  points  are  recorded. 


4.5  Conversion  to  Coefficient  Form 

The  aerodynamic  loads  acting  on  a  model  are  usually  expressed  in  the  form  of  coefficients, 
forces  being  represented  by 


and  moments  by 


(4.15) 


n 

qSb 


°  qSb'  Cm  qSc’  C" 

where  q  is  he  dynamic  pressure,  S  reference  area  and  h  and  c  are  reference  lengths. 


(4.16) 
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Dynamic  pressure  is  calc  'a ted  from  the  free  stream  static  pressure  \P)  and  the  free  stream 
total  pressure  ( H )  using  the  isotropic  flow  relationships 


giving 


and 


-  '} 


q  -  >,pl  -  -  1  yPM* 


With  y  —  I  '4  these  equations  become 


and 

q  --  0  7  PM-  (4.18) 

It  should  be  noted  that  the  assumption  (common  in  low  speed  wind  tunnels)  that 
q  H  P 

should  be  avoided  in  all  wind  tunnel  testing.  The  error  in  this  assumption  amounts  to  0-5",  at 
M  =  0- 14  (F  *  50  ms  ')  rising  to  more  than  2",,  at  M  0-3  <F  s  100  ms  '). 

For  many  tests  in  which  the  model  is  mounted  via  a  rear  sting,  the  pressure  on  the  blunt 
base  and  hence  the  axial  force  acting  on  the  model  is  influenced  to  a  large  extent  by  the  sting. 
For  tests  where  the  sting  enters  the  model  by  way  of  a  jet  engine  exhaust  duct,  this  problem 
may  be  overcome  by  correcting  the  axial  force  so  that  it  corresponds  to  the  zero  base  drag 
condition.  To  obtain  flight  values  of  axial  force,  the  engine  thrust  coefficient  is  then  added  to 
this  corrected  axial  force  to  give  the  total  axial  force.  For  models  without  jet  engine  exhausts, 
the  problem  of  obtaining  accurate  measurements  of  base  drag  is  much  more  difficult.  In  these 
cases  the  base  geometry  has  probably  been  grossly  distorted  to  allow  entry  of  the  sting.  The  base 
pressure  will  also  be  strongly  affected  by  the  presence  of  the  sting  and  by  the  condition  of  the 
model  boundary  layer.  The  only  practical  methods  for  estimating  base  drag  arc  to  do  a  separate 
series  of  tests  with  the  model  mounted  in  some  other  manner,  or  to  use  purely  empirical  methods. 
In  either  case,  results  from  sting  supported  tests  should  be  corrected  to  the  zero  base  drag  con¬ 
dition  to  provide  a  rational  starting  point. 

Hence  it  is  desirable  that  all  tests  of  sting  mounted  models  should  include  a  correction 
of  axial  force  to  the  zero  base  drag  condition.  Given  the  pressure  in  the  model  cavity  (fi,)  and 
the  area  of  the  base  "opening"  (Si,),  a  base  force  coefficient  may  be  defined  by 


Cxi, 


</’■,  P)Su 
qS 


(4.19) 


where  S  is  the  reference  area  used  to  non-dimensionalize  axial  force  coefficient.  Axial  force 
coefficient  corrected  to  zero  base  drag  i>  then  given  by 

Cx  Cx'  Cxh  (4.20) 
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where  Cx'  is  (he  axial  force  coefficit  .it  measured  by  the  balance. 

The  above  analysis  is  only  valid  for  cases  in  which  the  base  area  is  normal  to  the  balance 
axes  v-axis,  i.e.  perpendicular  to  .v n‘,  and  the  centroid  of  the  base  area  is  on  the  vB'  axis.  These 
conditions  are  often  satisfied  in  practical  models  (at  least  to  a  significant  level  of  approximation) 
However,  in  the  general  case,  base  force  will  contribute  to  loads  in  all  six  balance  components. 
For  this  case,  force  components  should  be  corrected  before  conversion  to  coefficient  form  by 


•Si.x 


F 


F'  (ffi,  P ) 


SnY 


(4.21) 


■S'tiz 


where  F"  are  the  uncorrccted  force  components  derived  from  balance  outputs  and  Stix,  Si>y 
and  Swi  are  the  projected  areas  of  the  base  in  the  ,vb  ,  .m  and  z»  directions  respectively. 

Similarly  the  three  moment  components  should  be  corrected  before  conversion  to  coefficient 
form  by 


0  c,  .iv 

•Si».x 

G  G  (/>,,  P) 

Zr  0  ,V<. 

•Shy 

.»•<•  x,  0 

St>z 

where  G  are  the  uncorrccted  moment  components  derived  from  balance  outputs  and  r,.. 
is  the  position  of  the  centroid  of  the  base  area  in  the  balance  axes  system. 

If.  as  may  occur  in  some  special  model  mounting  arrangements,  there  is  more  than  one 
base  area,  ihen  the  above  analysis  should  be  carried  out  for  each,  and  the  effects  summed. 

Where  hinge  moments  arc  measured  on  a  model,  it  is  customary  to  express  them  by  coefficients 
of  the  form 


C„ 


h 

<1  Sh  d h 


(4.23) 


where  h  is  the  measured  hinge  moment  and  5 h  and  d»  are  a  characteristic  area  and  length 
associated  with  the  control  surface. 


4.6  Transformation  to  Body  Axes 

Aerodynamic  coefficients  considered  thus  far  have  been  expressed  in  balance  axes.  The 
transformation  to  express  the  coefficients  in  body  axes  must  taVe  into  account  not  only  the 
possible  rotational  misalignment  of  the  two  axes  systems  but  also  the  non-coincidence  of  their 
origins. 

Consider  first  the  rotational  misalignment.  The  transformation  matrix  which  rotates  the 
balance  axes  into  the  body  axes  may  be  best  expressed  in  terms  of  three  readily  measured  angles 
defined  as : 

( i  -the  angle  between  the  projection  of  the  model  vu-axis  in  the  vbT  h  plane  and  the 
vn  -axis  positive  nose  to  starboard. 
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f-2 — the  angle  between  the  model  .VB-axis  and  its  projection  on  the  .vb)b  plane— positive 
nose  up. 

£3 — the  angle  between  the  model  zn-axis  and  its  projection  in  the  .vb'zb'  plane — positive 
clockwise  looking  forward. 

Note  that  this  set  of  angles  is  not  a  set  of  Euler  angles — each  rotation  is  about  the  undis¬ 
turbed  axis.  This  set  of  angles  is  chosen  however  to  make  measurement  of  the  required  angles 
as  simple  as  possible.  In  terms  of  these  angles  it  may  be  shown  that  the  transformation  matrix 
is  given  by 


[A]bb'  = 


cos  fi  COS  fi 

sin  sin  £3  —  A  cos  £1  sin  f>  cos-  fa 

ft 

cos  f i  sin  fi  cos2  fa  sin  fa  +  zf  sin  f  > 

B 


sin  f  1  cos  fa 


sin  f3 


—sin  f  > 

cos  f  >{eos  f  1  sin  £3^  A  sin  fi  sin  f-.( 

B 

cos  fal  sin  fi  sin  f •>  sin  (3  +  .4  cos  f  1  ;■ 

fi 


.  .  .  (4.24) 


where 


4  -  \  cos2  f3  sin2  fi  cos2  fa 

and 

B  =  cos2  f  1  cos2  fa  +  sin2  fa 

Hence  the  transformation  of  the  force  coefficients  is  given  by 

Fb  =  [A)bb  Fb’  (4.25) 

The  transformation  of  each  component  of  the  moment  must  also  include  the  contribution 
from  a  possible  shift  in  origin  between  the  two  axes  systems,  and  may  be  expressed  as 


0 


Gb  —  [AIbb*  Gb’  + 


0 


•Vli 


Fb 


(4,26) 


i'h 


b 


0 


where  (jri>,  >a,  zb)  is  the  location  of  the  origin  of  the  body  axes  system  with  respect  to  the  origin 
of  the  balance  axes  system,  b  is  the  reference  length  used  to  non-dtmensionalize  rolling-  and 
yawing-moment  coefficients,  and  <  the  reference  length  used  to  non-dimensionalize  the  pitching 
moment  coefficient.  Note  that  since  the  rotations  are  applied  first,  followed  by  the  translation. 
.xi>.  >'a.  zb  must  be  measured  in  body  (i.e.  model)  axes  rather  than  balance  axes. 
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5.  DETERMINATION  OF  MODEL  ATTITUDE 


Application  of  the  processes  described  in  the  previous  section  results  in  a  set  of  model 
aerodynamic  load  coefficients  expressed  in  the  body-axes  system.  Before  these  coefficients 
may  be  transferred  to  other  axes  systems,  the  aerodynamic  angles  *  and  /S  defining  the  direction 
of  the  relative  motion  between  the  model  and  the  wind  vector  must  be  determined.  This  in  turn 
requires  that  the  model  orientation  angles  p.  0  and  P  be  known. 

Practical  wind  tunnel  model  support  rigs  usually  only  allow  two  of  the  three  orientation 
angles  to  be  controlled.  The  two  most  popular  arrangements  are  the  yaw  pitch  {p.  0)  rig  and  the 
pitch  roll  t If  ii)  rig.  In  the  former,  the  model  is  yawed  through  an  angle  p  about  the  r»'  axis 
and  then  pitched  through  an  angle  t>  about  the  rotated  in'  axis.  In  the  latter  case,  the  model  is 
pitched  through  an  angle  »  about  the  i  n'  axis  and  then  rolled  through  an  angle  P  about  the 
rotated  .vif  axis.  In  both  cases  the  rotations  may  be  considered  as  sub-sets  of  the  general  rotations 
through  all  three  angles  {p.  «,  <i). 

In  the  past  several  years,  several  systems  have  been  proposed  for  measuring  the  orientation 
angles  directly."’111-  Most  of  these  systems  use  optical  sensors  of  some  type.  However,  until 
such  systems  are  in  general  use.  current,  less  satisfactory  approaches  must  be  used.  If  the  model 
is  large  enough,  the  preferred  approach  is  to  mount  three  accelerometers  inside  the  model  to 
measure  the  components  of  the  gravity  vector  in  three  orthogonal  directions.  If  the  model  is 
not  large  enough  to  accept  the  accelerometer  package  (approximately  75  mm  cube),  the  orien¬ 
tation  angles  must  be  computed  from  a  knowledge  of  the  attitude  of  the  sting  root  (i.e.  that  set 
by  the  model  support  rig).  In  this  case,  it  is  clear  that  any  deflection  of  the  balance  and  or  sting 
must  be  taken  into  account.  Although  not  as  obvious,  balance  and  or  sting  deflections  must 
also  he  taken  into  account  when  the  aerodynamic  angles  are  computed  from  the  outputs  of  a 
set  of  tri-axial  accelerometers. 

Finally,  in  obtaining  the  orientation  angles  (which  are  measured  relative  to  the  tunnel  axis 
system),  it  is  clear  that  any  angular  offsets  between  the  measuring  systems  in  use  and  the  tunnel 
axes  system  must  he  included. 


5.1  Calculation  of  Model  Deflections 

Differences  between  the  orientation  angles  and  the  angles  set  by  the  model  support  rig  will 
m  general  be  due  to  deflections  in  all  parts  of  the  mechanism  attaching  the  model  to  "earth". 
It  is  useful  to  divide  the  complete  mechanism  into  two  distinct  parts  when  considering  deflections. 
One  part,  including  the  balance  and  the  portion  of  sting  upstream  of  any  roll  angle  producing 
mechanism,  is  ammenable  to  a  general  analysis.  In  many  cases,  this  part  of  the  mechanism  can 
be  removed  and  its  deflection  characteristics  measured  experimentally.  The  remainder  of  the 
mechanism  cannot  be  easily  characterised  in  a  general  manner. 

For  the  removable  (balance  stingl  part  of  the  mechanism,  it  will  be  assumed  that  deflection 
angles  arc  small  enough  that  the  deflection  remains  elastic,  allowing  the  deflection  angles  to 
be  described  simply  by  a  linear  stiffness  times  a  load.  Hence  the  angular  deflections  may  be 
represented  as  functions  of  the  absolute  applied  loads  by 


|A)  Hn’ 


(5.1) 


where  u,  rj  and  x  are  the  Ismail)  angular  deflections  about  the  balance  axes  vn'.  iV  and  cb 
respectively  .  [A  ]  is  a  3  x6  matrix  of  stiffness,  and  Fin'  are  the  absolute  loads  applied  to  the  balance. 
Note  that  since  Hu'  are  absolute  loads,  including  contributions  From  both  aerodynamics  and 
gravity,  no  iteration  is  necessary  to  determine  the  angular  deflections  (as  would  be  the  case  if 
non-absolute  loads  were  used). 
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For  many  practical  balance/sting  arrangements,  many  of  the  36  stiffness  terms  in  [K]  will 
be  zero.  For  example,  if  the  balance/sting  neutral  axes  are  coincident  with  the  balance  axes 
(i.e.  a  fully  symmetric  arrangement),  then  [A"]  will  have  the  form 


0  0  0  *i  0  0 


[*]  = 


0 


0 


0 


0 


k,  0 


(5.2) 


j^O  0  0  0  0  a3J 

If,  as  is  common,  the  neutral  axes  coincide  with  the  xg'  and  i  n'  axes  but  not  the  ch'  axis,  two 
extra  terms  will  appear  in  [A]  thus 

fo  0  0  A,  0  0  1 


[*]  = 


0 


0 


a2  0 


A3  0 


(5.3) 


j^O  ki  0  0  0  AsJ 

The  deflection  angles  computed  from  (5.1)  will  be  expressed  in  the  balance-axes  system,  and 
must  be  transformed  to  the  body-axes  system  before  being  used  in  the  calculation  of  the  aero¬ 
dynamic  angles.  The  simplest  way  to  achieve  this  transformation  is  to  rewrite  (5.1)  as 


V  =  [*]Hb 


L.  J 


where  Hb  is  obtained  from 


(5.4) 


Hb  =  [3]bb'  Hb' 


(5.5) 


where  [A]Bb'  is  the  transformation  matrix  from  balance  to  body  axes  presented  in  (4.24). 

Deflection  of  the  part  of  the  model  support  mechanism  between  the  “earth"  end  of  the 
balance/sting  and  earth  depends  on  the  physical  layout,  and  hence  a  general  procedure  for  their 
treatment  is  not  possible.  In  most  if  not  all  cases,  this  part  of  the  support  will  have  been  designed 
to  be  much  stiffer  than  the  balance/sting  (for  which  the  physical  constraints  preclude  a  stiff 
design),  and  hence  its  contribution  to  the  overall  deflection  of  the  model  can  usually  be  ignored. 
Those  cases  for  which  this  assumption  is  not  valid  must  be  considered  on  an  individual  basis 
keeping  in  mind  that  some  parts  of  the  mechanism  may  move  while  others  remain  fixed  when 
setting  up  a  particular  model  attitude.  Deflections  from  all  parts  of  the  mechanism  must  be 
resolved  into  balance  axes  (and  then  to  body  axes)  for  each  attitude. 


5,2  Model  Attitude  from  Sting  Root  Attitude 

The  values  of  yaw,  pitch  and  roll  angles  indicated  by  the  model  support  rig  readouts  must 
be  corrected  for  angular  offsets  between  the  support  rig  zero  datums  and  the  tunnel  axes  system 


>!>  =  i 

e  =  eR+e0  >  (5.6) 

4  =  ^R  +  ^o  I 
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where  </ir,  Hr  and  i>&  are  the  values  indicated  by  the  support  rig  sensors  and  ifio.  do  and  4o  are 
angular  offsets  from  the  tunnel  axes  to  the  support  rig  zero  datums.  In  the  absence  of  significant 
deflections  in  the  non-removable  part  of  the  support  rig,  these  angles  represent  the  attitude  of 
the  sting  root.  With  the  addition  of  the  angular  deflections  of  the  balance  and  sting  (x,  rj,  v) 
computed  by  the  method  of  the  previous  section,  the  aerodynamic  angles  *  and  0  may  be  readily 
calculated.  Appendix  A  sets  out  the  required  formulae  ((A.  10)  and  (A.  12))  and  their  develop¬ 
ment.  It  should  be  noted  that  the  final  forms  of  these  formulae  presented  in  Appendix  A  are 
valid  only  while  the  angular  deflections  remain  small  (less  than  about  ±  2  if  the  desired  accuracy 
in  a  and  0  is  ±0  01  ).  This  will  be  the  case  for  almost  all  practical  sting/balance  combinations. 
Kor  other  cases,  a  more  detailed  analysis  will  be  required. 


5.3  Model  Attitude  from  Accelerometer  Outputs 

The  output  of  a  set  of  tri-axial  accelerometers  mounted  in  the  model  is  the  unit  gravity 
vector 


gA  - 


?XA 

{ya 

(z  A 


(5.7) 


whose  components  {xa.  {ya  and  £za  are  the  measured  components  of  the  unit  gravity  vector 
in  the  axes  system  of  the  accelerometer  mounting  block  <a  a.  .(  a,  r.x).  It  is  worth  noting  that  since 
gA  is  a  unit  vector  then 


v  Sxa'-Kya’+Sza2  =  '  (5.8) 

In  general,  the  axes  system  of  the  accelerometer  mounting  block  will  not  coincide  with  the 
body-axes  system.  Therefore,  to  be  useful  in  determining  the  model  attitude,  the  accelerometer 
outputs  must  first  be  transformed  to  the  body  axes  system  by 


gu  = 


Cxu 

Cyb 

{zb 


(3(iia  gA 


(5.9) 


The  transformation  matrix  [A]ba  which  rotates  the  accelerometer  mounting  block  axes 
into  the  body  axes  is  best  expressed  in  terms  of  three  readily  measured  angles  (which  are 
analogous  to  those  used  to  transform  the  balance  axes  to  body  axes)  defined  by 

on — the  angle  between  the  projection  of  the  model  vB-axis  in  the  va>  a  plane  and  the  ya 
axis — positive  to  starboard. 

un — the  angle  between  the  model  XH-axis  and  its  projection  on  the  ya.i  a  plane — positive 
nose  up. 

<«3 — the  angle  between  the  model  r«-axis  and  its  projection  in  the  jcaZa  plane— positive 
clockwise  looking  forward. 
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The  transformation  matrix  [A)ha  is  then  given  by  equation  (4.23)  with  (fe,  {•>,  fe)  replaced 
by  iwi,  wj,  wa).  Again  it  should  be  noted  that  (on.  u>>,  0.3)  is  not  a  set  of  Euler  angles.  The  angles 
(o>i.  tu2,  013)  can  be  simply  determined  by  measuring  the  attitude  of  the  model  body-axes  system 
with  respect  to  the  tunnel-axes  system  when,  with  the  wind  otf,  the  output  of  the  accelerometers 
is  given  by 


fetA 

0 

fe'A 

0 

fee 

1 

_  _ 

(5.10) 


i.e.  when  the  accelerometer  axes  are  aligned  with  the  gravity  axes. 

The  pitch  and  rod  angles  of  the  body-axes  system  will  then  be  given  by 


arcsin  (ixn)  +  H„ 

(5.11) 

arCtan(fe») 

(5.12) 

where  fl0  is  the  angular  offset  between  the  tunnel  free  stream  velocity  vector  and  the  horizontal 
plane. 

For  the  large  pitch  angles  (H  >  v  4)  it  will  be  more  accurate  to  determine  the  pitch  angle 

from 


I*  •-  arccos  x  Sv it2+fe.ir+  (5.13) 

It  should  also  be  noted  that  as  the  pitch  angle  appn  tches  n;2.  the  feiiand  fen  axes  approach 
the  horizontal  plane  and  changes  in  roll  angle  will  no  longer  produce  changes  in  the  orientations 
of  these  accelerometers  with  respect  to  the  gravity  vector.  Hence  as  the  pitch  angle  approaches 
jt/2.  all  roll  angle  information  is  inevitably  lost,  and  the  accuracy  of  the  roll  angles  derived  from 
(5.12)  as  the  pitch  angle  approaches  n  2  will  be  significantly  degraded.  For  a  particular  test 
where  these  limitations  are  important,  it  is  of  course  possible  to  change  the  orientation  of  the 
accelerometer  package,  or  to  add  extra  accelerometers  to  ensure  an  accurate  determination  of 
roll  angle. 

Although  these  values  of  pitch  and  roll  angles  include  the  effects  of  the  deflection  of  the 
support  rig.  the  accelerometer  package  is  insensitive  to  rotations  in  the  horizontal  plane,  and 
hence  to  any  such  rotation  produced  by  support  rig  deflection.  This  rotation,  which  for  non¬ 
zero  roll  angles  will  include  contributions  from  angular  deflections  about  both  the  rotated  zn-axis 
and  the  rotated  in-ax  is,  may  be  evaluated  to  sufficient  accuracy  (so  long  as  angular  deflections 
remain  small)  from  the  angular  deflections  calculated  by  the  methods  of  section  5.1  by 


tan  if* 


sin  (x  cos  if  ysin*) 
cos  »  cos  (x  cos  </>  r/  sin  if) 


(5  14) 


where  l>  and  if  are  those  determined  from  the  accelerometer  package.  The  yaw  angle  of  the  body- 
axes  system  will  then  be  given  by 


if  ifK+if*  +  0„ 


(5.15) 


where  if  r  is  the  yaw  angle  indicated  by  the  support  rig  sensor  and  1 f„  is  the  angular  offset  between 
the  support  rig  yaw  angle  datum  and  the  v-r-axis. 

The  aerodynamic  angles  *  and  may  now  be  calculated  from  if.  »  and  if  using  equation  ( A .  I ). 
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6.  CORRECTIONS  FOR  THE  EFFECT  OF  WIND-TUNNEL  WALL  INTERFERENCE 


The  measurement  of  vehicle  characteristics  from  wind  tunnel  tests,  especially  at  transonic 
speeds,  is  subject  to  several  sources  of  experimental  error.  One  of  the  most  significant  of  these 
is  wind  tunnel  wall  interference.  This  section  addresses  the  problem  of  correcting  experimental 
data  for  the  effects  of  wind  tunnel  wall  interference. 

When  a  model  is  placed  in  a  wind  tunnel,  the  aerodynamic  loads  measured  on  it  differ  from 
those  which  would  be  measured  if  the  model  were  tested  under  similar  conditions  in  free  air, 
because  the  wind  tunnel  walls  interfere  with  the  flow  around  the  model.  It  is  usual  to  make  correc¬ 
tions  to  measured  quantities,  and  to  Mach  number  and  model  attitude,  so  that  the  corrected 
results  correspond  to  those  which  would  be  obtained  if  the  model  were  tested  in  free  air  at  the 
corrected  Mach  number  and  attitude.  The  magnitude  of  these  corrections  may  be  reduced  to 
arbitrarily  small  values  by  using  models  which  are  small  when  compared  to  the  test  section. 
However,  in  an  effort  to  make  maximum  possible  use  of  the  available  unit  Reynolds  number, 
and  hence  to  minimize  errors  due  to  scale  effects,  models  are  generally  designed  with  a  size  which 
requires  significant  corrections  for  the  effects  of  wind  tunnel  wall  interference. 

For  a  model  of  reasonable  size  (compared  to  the  test  section),  these  corrections  may  be 
calculated  from  linearized  potential  flow  theory.13  Changes  in  local  stream  direction  and  stream 
curvature  caused  by  the  presence  of  the  wind  tunnel  walls  preventing  the  normal  flow  curvature 
over  a  lifting  body  that  occurs  in  free  air,  are  referred  to  as  lift  interference.  Changes  in  local 
stream  velocity  and  longitudinal  gradients  caused  by  the  lateral  constraint  of  the  flow  past  the 
model  are  referred  to  as  blockage  interference.  The  two  types  of  interference  are  assumed  to  be 
mutually  independent. 

This  classical  approach  to  estimating  the  effects  of  wind  tunnel  wall  interference  suffers 
from  inaccuracy  due  to  its  first  order  linearized  approach  and  the  necessarily  simple  way  in  which 
the  physical  model  is  mathematically  represented.  Despite  these  limitations,  for  tests  in  tunnels 
with  solid  walls,  the  results  produced  are  sufficiently  accurate  for  most  development  testing. 
The  results  for  tests  in  transonic  tunnels,  which  are  usually  equipped  with  ventilated  (either 
porous  or  slotted)  walls,  are  however  far  from  satisfactory.  The  problem  here  is  one  of  correctly 
representing  the  boundary  conditions  at  the  tunnel  walls.  In  the  case  of  solid  walls,  the  boundary 
conditions  are  straightforward  and  mathematically  amenable.  For  ventilated  walls,  the  wall 
boundary  conditions  are  complex,  and  it  is  usual  to  change  the  mixed  (open,  closed)  boundary 
conditions  into  an  equivalent  homogeneious  one,  which  produces  the  same  solution  at  the  centre 
of  the  test  section,  but  the  correct  representation  remains  open  to  question. u’15  In  most  cases, 
recourse  must  be  made  to  experiment  to  determine  the  correct  wall  boundary  conditions  for  each 
particular  test  section. 

In  recent  years  it  has  been  realized  that  the  classical  homogenious  boundary  conditions  for 
ventilated  walls  is  not  satisfactory,  and  further  that  classical  linear  theory  itself  is  not  adequate 
to  describe  the  truly  non-linear  nature  of  flows  at  supercritical  speeds.  In  view  of  the  fact  that 
the  boundary  conditions  at  the  wall  are  not  well  defined  and  are  in  fact  model  dependent  (through 
the  pressure  and  flow  fields  imposed  at  the  wall),  the  use  of  measured  wall  pressure  data  to  infer 
the  wall  constraints  would  appear  to  be  a  logical  approach.  Many  proposals  to  make  use  of  wall 
pressure  data  in  computing  interference  corrections  have  therefore  been  put  forward.  The  most 
promising  of  these  also  take  advantage  of  advances  in  our  ability  to  compute  transonic  flows, 
to  combine  a  truly  non-linear  calculation  method  with  measured  wall  pressure  distributions. 
Such  methods  arc  now  fairly  well  developed  for  the  two-dimensional  aerofoil  case,16’17  but  their 
routine  application  to  the  testing  of  three-dimensional  models  is  still  some  way  off.18  Until 
that  time,  the  estimation  of  wall  interference  effects  in  ventilated  test  sections  must  rely  on  the 
questionable  methods  of  classical  theory. 


6.1  Transformation  to  Wind  Axes 

Classical  linear  interference  theory  has  been  developed  to  make  use  of  measured  forces  and 
moments  expressed  in  the  wind-axes  system.  Hence  before  applying  any  interference  corrections, 
forces  and  moments  must  be  transformed  from  the  body-axes  system  to  the  wind-axes  system 


Fw  —  [A)wb  Fb 


(6.1) 


and 


Gw  —  [A]wb  Gb  (6.2) 

where  [A]w  B  is  the  required  transformation  matrix.  The  required  matrix  is  presented  in  Appendix 
B  in  terms  of  all  common  combinations  of  the  aerodynamic  angles  x..  (3»,  on  and  fit- 


6.2  Application  of  Classical  Theory  is  Solid-Wall  Tunnels 

For  tunnels  with  solid  walls,  the  wall  boundary  conditions  are  well  defined  and.  at  least  for 
subcritical  flows,  the  results  of  classical  interference  theory  may  be  used  for  correcting  experi¬ 
mental  measurements. 

Consider  first  the  corrections  which  must  be  made  for  the  lateral  constraint  of  the  flow 
past  the  model  due  to  the  presence  of  the  walls,  i.e.  the  blockage  corrections.  This  effect  is 
equivalent  to  an  increase  in  the  free  stream  velocity  in  the  vicinity  of  the  model,  and  for  models 
of  normal  size,  may  be  considered  constant  over  the  model.  Hence  the  corrected  free  stream 
velocity  is  given  by 


ft  f't(l  +  f)  (6.3) 

where  <,  the  blockage  factor,  is  usually  split  into  separate  contributions  from  the  body  and  from 
its  wake 


:  (6.4) 

For  a  three-dimensional  aircraft  model,  the  contribution  from  the  body,  the  solid  body 
blockage,  <„  can  be  further  split  into  contributions  from  the  wing  and  from  the  fuselage 


r~  *nt)  T  r*w 


(6,5) 


According  to  classical  theory,  the  contribution  from  the  wing  may  be  represented  by 

where  C  is  the  tunnel  cross-sectional  area.  F»  is  the  volume  of  the  wing  and  T  is  a  constant  which 
depends  on  the  tunnel  cross-sectional  shape  ( B!H )  and  the  ratio  of  wing  span  to  tunnel  width 
(, h/B ).  For  a  streamlined  body  of  revolution,  with  length  to  diameter  ratio  /  the  solid  body 
blockage  is  given  by 


T  I  /  0-4/3\ 

e»-  =  c3ia^('+  f 

where  Ft,  is  the  volume  of  the  body. 

The  contribution  of  the  wake  to  the  blockage  factor  is  given  by 


(6.7) 


— 


S 

4c 


I  +0  4A/,2 

P 


Cv0 


(6.8) 


where  S  is  the  wing  area  and  M,  the  (uncorrected)  free  stream  Mach  number,  and  Cd0  is  the  model 
zero-lift  drag  coefficient. 

With  contributions  to  the  blockage  factor  c  calculated  from  (6.6),  (6.7)  and  (6.8),  free  stream 
quantities  may  now  be  corrected  by  (for  y  =  14): 
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qt  —  {1  +(2  -  Mchiqi 

(6.9) 

|)f=jl  M,2f)p, 

(6.10) 

Mi  =  {1  +(l  +0  2W,2)e}W, 

(6.11) 

Ret  —  { 1  +( 1  01  MrU}Ret 

(6.12) 

Note  that  classical  thoery  makes  no  allowance  for  bluff  body  (i.e.  stalled  wing)  blockage. 
This  contribution  will  be  significant  in  the  post-stall  region  where  it  may  produce  an  additional 
correction  to  dynamic  pressure  on  the  order  of  I 

Turning  now  to  lift  interference.  The  presence  of  the  upper  and  lower  wind  tunnel  walls 
modifies  the  flow  curvature  over  a  lifting  body  that  occurs  in  free  air.  Corrections  must  therefore 
be  made  to  the  angle  of  incidence  and  to  measured  coefficients  for  the  mean  value  and  streamwise 
variation  of  the  induced  upwash.  By  choosing  the  location  of  the  vortex  used  by  classical  theory 
to  represent  the  circulation  to  be  at  the  centre  of  the  wing  (i.e.  mid  chord),  all  the  upwash  correc¬ 
tions  may  be  included  in  corrections  to  the  angle  of  incidence,  there  being  no  correction  needed 
for  lift  coefficient.  There  remains,  however,  a  residual  correction  to  pitching  moment  coefficient. 
The  correction  to  angle  of  incidence  will  in  this  case  be  given  by 


-- 


*<  + 


S 
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So  + 


Ac 

2fJ  H 


*) 


(6.13) 


where  8,,  is  related  to  the  mean  value  of  the  induced  upwash,  and  &i  to  the  streamwise  gradient 
i.e.  to  streamline  curvature.  Note  the  relative  increase  in  importance  of  the  streamline  curvature 
term  with  increasing  Mach  number  (decreasing  (3).  It  has  been  usual  (especially  in  low-speed 
tunnel  work)  to  express  equation  (6.13)  in  the  form 


if  --  +■  S„ 


s 

c 


(6.14) 


8„  is  then  a  function  only  of  tunnel  shape  and  may  be  determined  either  experimentally  or 
theoretically  for  a  particular  tunnel.  The  constant  r-<  on  the  other  hand  is  a  function  of  both 
model  and  tunnel  geometry  (and  Mach  number). 

The  residual  correction  to  pitching  moment  coefficient  is  given  by 


Cmf 


r  U*)s< 

Vnil  +  ^  (1 


(6.15) 


where  a  is  the  wing  lift-curve  slope,  and  (Aa)„,.  is  the  streamline  curvature  induced  correction 
to  angle  of  incidence  which,  from  equation  (6.13)  is  given  by 


(•ia).,. 


S 

lVHh'cC'' 


or  from  equation  (6.14) 


0*1-, 


(6.16) 


(6.17) 


The  upwash  induced  by  the  tunnel  walls  also  has  the  effect  of  reducing  the  measured  drag 
due  to  the  wing  appearing  to  have  a  higher  aspect  ratio.  According  to  classical  theory  the  correc¬ 
tion  to  drag  coefficient  is  given  by 


Cor 


On  +  So 


s 

c 


Cl.,2 


(6.18) 
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Similarly  the  induced  upwash  will  influence  the  flow  field  of  any  tailplane  which  must  also 
be  corrected.  The  correction  to  the  tailplane  angle  takes  the  form  of  equations  (6.13)  or  (6.14) 
with  C/2  replaced  by  It  (the  distance  of  the  tailplane  aerodynamic  centre  aft  of  the  aerodynamic 
centre  of  the  wing)  in  (6.13)  and  with  r2  evaluated  at  the  tailplane  location  in  (6.14). 

An  additional  correction  is  required  to  the  drag  coefficient  because  the  images  of  the  source 
representing  the  wake  impose  a  longitudinal  pressure  gradient  along  the  tunnel  and  hence  a 
longitudinal  bouyancy  force  on  the  model.  This  pressure  gradient  may  conveniently  be  regarded 
as  linear  along  the  length  of  the  model.  This  correction  will  then  be  given  by 

ACd  =  C(k<!+0  4A/,2)(„  (6.19) 

The  total  correction  to  drag  coefficient  for  the  effects  of  both  lift  and  blockage  interference 
and  for  longitudinal  bouyancy  therefore  becomes 

Cor  =  jcnt  +  8„  £  CL,2  +  (I  +  0'4M,2)CDlj/jl  +(2-  Mr’)  (6.20) 

and  the  total  correction  to  pitching  moment  coefficient  becomes 

-r}/{,+ 

For  all  other  coefficients,  the  only  correction  required  is  that  due  to  blockage  so  that  lift 
coefficient,  for  example,  becomes 


- 


Cmf  — 


(2~M,2)  €} 


(6.21) 


Cu  =  Cu/{l  +(2- M\2)e]  (6.22) 

The  results  from  classical  interference  theory  quoted  above  are  not  valid  in  the  most  general 
case.  Several  simplifying  assumptions  have  been  made  as  follows : 

(i)  It  has  been  assumed  that  the  model  and  specifically  the  centre  of  pressure  of  the  wing 
remain  on  the  tunnel  centre  line.  Movements  off  the  centre  line  will  require  modifi¬ 
cations  to  both  the  lift  interference  factors  So  and  Si  and  the  solid  blockage  factor 

(ii)  No  account  has  been  taken  of  effects  due  to  asymmetric  wing  loading  with  the  model 
symmetric  in  the  tunnel  nor  to  asymmetric  placement  of  the  model  in  the  tunnel  (both 
due  to  rolling  and  yawing  of  the  model). 

(iii)  No  account  has  been  taken  of  possible  bouyancy  effects  (affecting  drag  coefficient) 
due  to  a  variation  of  the  empty  tunnel  longitudinal  pressure  along  the  test  section  axis. 


In  general  these  assumptions  produce  small  errors,  within  the  expected  accuracy  of  classical 
theory.  All  of  these  effects  can  however  be  included  in  the  structure  of  the  theory,  and  if  any  of 
the  above  assumptions  are  severely  violated,  the  magnitude  of  the  error  should  be  investigated 
and  included  if  significant. 


6.3  Application  of  Classical  Theory  in  Ventilated-Wall  Tunnels 

In  principle,  the  walls  of  a  ventilated  wall  tunnel  affect  the  flow  over  a  model  in  a  similar 
manner  to  those  of  a  tunnel  with  open  or  closed  walls.  Since  many  interference  parameters 
change  sign  between  the  open  and  closed  wall  cases,  it  is  theoretically  possible  to  design  ventilated 
walls  which  reduce  or  eliminate  some  sources  of  interference.  However,  in  practice  ventilated 
walls  are  often  arranged  primarily  for  the  generation  of  low  supersonic  Mach  numbers  using 


f  ; — 
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diffuser  suction,  or  to  minimize  shock  reflections  at  high  subsonic  and  low  supersonic  Mach 
numbers,  with  little  attempt  made  to  utilize  the  available  cancellation  of  interference. 

The  application  of  classical  interference  theory  to  wind  tunnels  with  ventilated  walls  makes 
use  of  the  same  set  of  relationships  as  presented  for  solid-wall  tunnels  in  the  previous  section, 
but  with  differing  values  of  the  parameters  «,  K  and  Si.  The  approach  of  linear  interference  theory 
to  calculating  these  parameters  is  to  represent  the  ventilated  walls  by  an  equivalent  homogenious 
boundary  conditions  characterised  (for  slotted  walls)  by  two  wall  parameters.  The  first  of  these, 
the  slot  parameter  F,  is  a  function  of  the  geometry  of  the  slotted  wall  given  by 

2d  (  Tra) 

fvn  {«■*  4  <6n> 

where  d  is  the  periodic  spacing  of  the  slots  whose  width  is  a.  F  is  supposed  to  represent  the  effect 
of  the  slot  geometry  on  the  tunnel  centre  line.  Equation  (6.23)  is  a  very  simple  representation, 
taking  no  account  of  the  finite  thickness  of  the  slots,  or  the  separation  on  the  plenum  side  of  the 
slot.  Many  more  complex  expressions  (see  for  example  Refs  14  and  15)  have  been  proposed  to 
better  represent  the  slot  geometry. 

The  second  parameter  characterising  the  slotted  wall,  the  porosity  parameter  P  (or  more 
usually  ft  P).  is  also  used  to  characterise  the  porous  wall.  It  is  introduced  in  the  slotted  wall  case 
to  take  account  of  viscous  effects  in  the  slots,  and  is  analogous  to  the  term  expressing  the  linearized 
pressure  drop  across  a  porous  wall.  Unlike  F.  there  is  no  way  of  computing  ft  P  a  priori,  but 
the  effects  of  viscosity  on  the  flow  through  a  slot  are  obviously  important,  especially  when  the 
slot  width  is  on  the  same  order  as  the  wall  boundary  layer  thickness.  It  is  doubtful  however, 
even  if  ft  P  correctly  characterizes  the  outflow  of  working  section  fluid,  whether  a  simple  linear 
parameter  can  ever  hope  to  characterize  the  inflow  of  fluid  from  the  plenum  chamber.  In  this 
case  it  is  possible  that  the  fluid  entering  the  slots  could  be  more  or  less  stagnant.  The  situation 
is  further  complicated  by  the  influence  of  the  model  on  the  growth  of  tunnel  wall  boundary 
layers,  in  w  hich  case  it  becomes  doubtful  whether  a  single  value  of  ft/P  could  be  valid  for  all 
tests  with  particular  walls,  or  even  the  whole  of  one  wall  for  a  particular  test. 

Given  these  problems,  it  is  not  surprising  that  classical  linear  theory  is  much  less  successful 
in  predicting  interference  in  ventilated  wall  tunnels  than  in  solid  wall  tunnels.  Even  for  the 
mathematically  simpler  case  of  the  two-dimensional  aerofoil,  experiments19  confirm  the  basic- 
failure  of  linear  theory  to  correctly  predict  interference  effects.  For  three-dimensional  testing 
the  results  arc  even  less  satisfactory.  Hence  it  has  been  the  practice  at  A  R  E.  not  to  correct 
results  obtained  in  ventilated  (slotted)  walls  for  interference.  This  approach  is  justified  on  the 
basis  of  the  following: 

(i)  Ventilated  walls  are  used  only  for  testing  at  high  (transonic)  speeds  for  which  other 
considerations  encourage  the  use  of  models  which  are  smaller  (when  compared  to  test 
section  area)  than  are  generally  used  at  low  speeds,  reducing  the  magnitude  of  the  inter¬ 
ference  effects. 

(ii)  The  ventilated  walls  in  use  in  A.R.L.  tunnels,  although  not  designed  to  eliminate 
interference,  do  produce  significantly  smaller  interference  than  do  solid  walls. 

(iii)  The  probable  inaccuracies  in  linear  interference  theory  for  ventilated  walls  are  expected 
to  be  on  the  same  order  as  the  total  interference  corrections 


6.4  Transformation  to  Other  Axes  Systems 

It  is  not  unusual  for  wind  tunnel  test  customers  to  request  that  data  be  presented  in  an  axes 
system  other  than  the  wind-axes  system.  Hence  once  interference  corrections  have  been  applied 
the  data  must  be  transformed  back  to  the  body-axes  and  or  stability-axes  system.  This  trans¬ 
formation  is  achieved  via  the  application  of  a  suitable  transformation  matrix,  for  example  from 
wind-axes  to  stability-axes  via 
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F  =  [A]sw  F w 


and 


Gs  =  (Ajsw  Gw 

Transformation  matrices  to  transform  data  between  body,  stability  and  wind  axes  are  pre¬ 
sented  in  Appendix  B  in  terms  of  all  common  combinations  of  the  aerodynamic  angles  *».  ft, 
»t  and  ft 


7.  CONCLUSIONS 

This  report  has  presented  a  general  approach  to  the  reduction  and  processing  of  wind 
tunnel  data.  It  has  not  attempted  to  provide  a  general  approach  for  all  possible  types  of  wind 
tunnel  testing,  but  has  concentrated  on  procedures  applicable  to  force  and  moment  data  acquired 
from  internal  sting  mounted  strain  gauge  balances.  The  procedures  should  however  be  applicable, 
with  minor  variations,  to  force  and  moment  data  acquired  from  modern  load-cell  based  external 
balances.  Although  every  attemp  has  been  made  to  keep  the  approach  as  general  as  possible 
within  its  area  of  application,  not  all  possibilities  have  been  included.  For  example,  the  practice, 
found  necessary  at  some  wind  tunnels,  of  using  separate  strain  gauge  balance  calibration  factors 
for  each  sign  of  applied  load,  has  not  been  considered.  This  practice  has  so  far  been  found  to  be 
unnecessary  with  the  types  of  balance  in  use  at  A.R.L.  Such  decisions  have  also  been  made  in 
several  other  areas.  Overall,  the  procedures  as  presented  are  sufficiently  general  to  include  all 
those  effects  which  have  been  found  necessary  for  the  types  of  tests  and  balances  so  far  encountered 
in  A.R.L.  wind  tunnels. 

Finally,  it  is  worth  noting  once  again,  that  no  matter  how  general  (and  hopefully  accurate) 
the  computer  representations  of  the  algorithms  included  in  a  data  acquisition  system  may  be. 
the  final  results  will  only  be  as  accurate  as  the  least  accurate  information  used  in  the  reduction 
process.  Thus,  no  matter  how  good  the  data  acquisition  system,  an  inaccurate  balance  cali¬ 
bration  matrix  for  example,  will  ensure  inaccurate  results.  The  often  quoted  adage  concerning 
computer  processing  —  garbage  in.  garbage  out  —  has  never  been  more  appropriate. 
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APPENDIX  A 


Relationships  Between  Aerodynamic  and  Orientation  Angles 

When  a  wind  tunnel  model  is  set  up  at  some  attitude,  the  set  of  Euler  angles  (</>,  0, p)  through 
which  the  model  is  rotated  to  reach  that  altitude  may  be  related  to  the  aerodynamic  angles  x 
and  £  by 

sin  as  =  cos  p  sin  0  cos  p  +  sin  p  sin  p 
sin  ft  =  cos  p  sin  0  sin  <p  -sin  ip  cos  <p 


tan  w>  =  tan  0  cos  <p  + 


tan  ip  sin  <p 


tan  ft  =  tan  0  sin  p  — 


tan  ip  cos  < 


It  should  be  noted  that  the  order  of  the  rotations  ( ip ,  0,  <p)  is  significant.  If  another  of  the 
six  possible  rotation  orders  is  used,  then  the  above  relationships  no  longer  hold.  For  example, 
if  the  rotations  are  applied  in  the  order  (0,  <p,  ip),  then  the  above  relationships  become 

sin  on  =  sin  ft  cos  p 

sin  ft  =  sin  ft  sin  p  cos  ip  -cos  ft  sin  ip 

tan  0  cos  <P  (A  2) 

tan  on  , 

cos  ip  +  tan  0  sin  <p  sin  i P 


tan  ft  = 


tan  0  sin  <£-tan  ip 
I  +  tan  ft  sin  <p  tan  ip 


In  most  wind  tunnels,  model  support  mechanisms  are  capable  of  providing  only  a  subset 
of  the  full  ( ip ,  0,  <P)  set  of  Euler  rotations.  The  most  commonly  available  rigs  are  the  pitch/roll 
(0,  <p)  rig,  and  the  yaw/pitch  rig  (p,  ft).  Such  rigs  produce  the  subset  of  Euler  rotations  in  the  same 
order  as  the  full  set.  In  such  cases  the  general  relationships  (A.I)  simplify  considerably.  For  the 
pitch/roll  (0,  P)  rig  they  become 

sin  oc„  =  sin  ft  cos  <p 
sin  ft  -  sin  0  sin  <p 

*  (A. 3) 

tan  *i  =  tan  0  cos  P 
tan  ft  =  tan  0  sin  p 


and  for  the  yaw /pitch  [p,  ft)  rig 
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sin  a6  =  sin  0  cos  4> 
sin  ft  =  —sin  0 
tan  «t  =  tan  9 
tan  ft  =  —tan  0/ cos  0 


(A.4) 


The  last  set  explain  the  reason  for  the  popularity  of  at  and  ft  as  the  choice  of  aerodynamic 
angles  in  low-speed  wind  tunnels  (where  model  rigs  are  traditionally  of  the  yaw/pitch  type), 
for  in  this  case 


ft  =  -0 

aj  —  9 


(A. 5) 


As  well  as  allowing  the  use  of  the  above  simplified  relationships  the  limitation  to  two 
Euler  rotations  also  makes  available  some  particularly  simple  inverse  relationships.  These  are, 
for  the  pitch/roll  (9  </>)  rig 


cos  9  =  cos  at  cos  ft 
tan  0  —  tan  ft/sin  a, 


sin2  0  =  sin2  a»  +  sin2  ft 
tan  0  —  sin  ft/sin  «s 

tan2  0  =  tan2  at  +  tan2  ft 
tan  0  =  tan  ft/tan  at 


(A.6) 


and  for  the  yaw/pitch  (0  0)  rig 


9  =  at 


sin  0  -=  sin  a,/cos  ft  . 

,  „  r-  <A-7> 

9  =  -ft 

9  —  at 

tan  i/i  =  —tan  ft/cos  at 


Note  that  the  above  inverse  relationships  are  grouped  into  three  pairs  of  equations  for  9 
and  ^  or  0  and  9,  one  pair  for  each  possible  choice  of  pairs  of  aerodynamic  angles  —  (at,  ft), 
(*•.  ft)  or  («t,  ft).  It  is  apparent  that  for  the  pitch/roll  rig,  the  above  expressions  introduce  an 
ambiguity  in  the  correct  sign  of  0.  However,  if  the  roll  angle  is  in  the  range  —  w/2<^jt/2 
(a  range  which,  when  combined  with  pitch  angles  in  the  range  —n/2<9siir/2,  gives  access  to 
all  values  of  a  and  ft  then  the  correct  sign  of  9  is  the  same  as  that  of  a  (either  a„  or  a().  Similar 
inverse  relationships  cannot  be  derived  for  the  general  (0,  0,  0)  case  since  in  that  case  any  given 
combination  of  aerodynamic  angles  can  be  produced  by  an  infinite  number  of  combinations 
of  the  three  orientation  angles. 
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The  relationships  presented  thus  far  should  be  interpreted  as  relating  the  aerodynamic 
angles  and  the  orientation  angles  of  the  sting  root,  that  is  the  orientation  angles  as  set  by  the 
model  support  rig.  In  practice  the  sting  and  balance  will  deflect  under  the  applied  loads  (both 
gravitational  and  aerodynamic)  and  the  orientation  angles  of  the  body-axes  system  will  differ 
from  those  of  the  sting  root.  Unless  model  attitude  is  measured  via  on-board  accelerometers, 
these  deflections  must  be  taken  into  account  when  calculating  the  aerodynamic  angles.  (See 
Section  5  for  the  effect  of  deflections  on  accelerometer  measured  attitudes).  The  effect  of  de¬ 
flections  on  the  relationship  between  aerodynamic  and  orientation  angles  for  both  the  pitch/roll 
rig  and  the  yaw/pitch  rig  are  presented  below.  In  the  derivation  of  these  relationships  it  has  been 
assumed  that  the  deflections  are  small,  and  hence  that  the  order  of  application  of  the  rotations 
is  not  significant.  Hence  the  order  of  application  chosen  for  each  case  is  that  which  results 
in  the  simplest  relationships.  However,  it  is  worth  noting  that  the  relationships  are  exact  for  the 
set  of  Euler  rotations  applied  in  the  stated  order. 

Consider  first  a  pitch/roll  rig  (ft  &)  for  which  the  deflections  of  the  balance  and  sting  are 
given  by 

X — deflection  in  the  yaw  direction  —  i.e.  a  rotation  about  the  0 z  axis. 
rj — deflection  in  the  pitch  direction  —  i.e.  a  rotation  about  the  0.y  axis. 
v — deflection  in  the  roll  direction  —  i.e.  a  rotation  about  the  Ox  axis. 


Then  for  deflections  applied  in  the  order  ( u.  *),  the  aerodynamic  angles  are  given  by 


sin  «s  =  sin  0  cos  <f>'  cos  ij  +  cos  8  sin  n 

sin  ft  =  sin  6  sin  4'  cos  x  -sin  x  (cos  @  cos  T  — sin  8  sin  i;) 

sin  8  cos  4‘  cos  n  +  cos  8  sin  n 
fan  *<  =  .  „  '■■■'.  ,,  . 

cos  *  (cos  0  cos  sin  8  cos  4  sin^  +  sinfisin^  sin  x  - 


(A.8) 


sin  0sin  4'  cos  x~ sin  x  (cos  0  cos  »?-sin  8  sin  i?) 
cosx(cos  0 cos y  —  sin  dcos^'sin  >;)  +  sin  0sin<ft  sin  x 


where 


4'  =  <j>+  u 


In  most  cases,  the  deflections  are  small  enough  that  the  approximations 
cos  x  cos  >;  =  cos  v  %  I 
sin  x  *  x !  sin  »  ip ;  sin  v  «  u 
can  be  made  and  (A.8)  can  be  reduced  to 

sin  =  sin  0  cos  4'  +  tj  cos  8 
sin  ft  =  sin  8  sin  4'  —  X  cos  0 


tan  0cos  4'  +  v 
I  +tan  8  (x  sin  4'  +  tj  cos  </>') 

tan  Ssin  4'  ~x 
1  +  tan  8  (x  sin  4'  -t)  cos  4 ') 

4'  —  4  +  u 


tan  = 


tan  ft  = 


(A. 9) 


(A. 10) 


f 


For  a  yaw/pitch  (0,  8)  rig,  applying  the  same  deflections,  but  in  the  order  (ij,  v,  41  gives 


APPENDIX  B 


Axes  Transformations 

This  Appendix  presents  a  collection  of  transformation  matrices  for  the  transformation 
of  vectors  between  the  body-axes  system  and  both  the  stability-  and  wind-axes  systems.  Trans¬ 
formations  are  given  in  terms  of  the  matrix  [A]i2  where  the  components  of  a  vector  A  in  axes- 
system  I  are  rotated  to  those  in  axes-system  2  by 


or 


Axi 

Axi 

Art 

=  [A]l2 

A  Y2 

Azi 

Azz 

Ai  =  [A]i2  A2 


The  matrix  [A]i2  is  given  in  terms  of  all  the  useful  combinations  of  the  aerodynamic  angles 
«s,  ft,  <*t  and  ft. 


B1  .—Transformations  in  Terms  of  «».  ft- 
From  body  axes  to  stability  axes. 


[A]sb  = 


cos  at  0  sin  at 

0  1  0 

—  sin  a(  0  cos  at 


From  body  axes  to  wind  axes : 


(A]wb  — 


cos  a(  cos  ft  sin  ft 
-cos  a(  sin  ft  cos  ft 
-  sin  a(  0 


"VT  - 


sin  at  cos  ft 
-sin  at  sin  ft 
cos  at 


B2. — Transformations  in  Terms  of  a8,  ft. 


From  body  axes  to  stability  axes : 


sec  ft,  \/cos2  aa  —  sin2  ft  0 


[■Msb  — 


0 


I 


—sin  as  sec  ft 


0 


sin  a8  sec  ft 
0 

sec  ft  v'cos2  ota  —  sin2  ft 


From  body  axes  to  wind  axes : 


Mwb  = 


V  cos2  as  -  sin2  ft 


tan  ft  v  cos2  ots  —sin2  ft 
—sin  aa  sec  ft 


sin  ft  sin  as 

cos  ft  —sin  ot8  tan  ft 

0  sec  ft\/cos2  otg  —  sin2  ft 


B3. — Transformations  in  Terms  of  at,  ft. 
From  body  axes  to  stability  axes : 


Msb  = 


cos  at 
0 

-sin  at 


sin  at 
0 

COS  at 


From  body  axes  to  wind  axes: 


I'Mw'B 


1 

s  sec2  at  +  tan2  ft 

cos  »t  tan  ft 
x  sec2  at  +  tan2  ft 

sin  at 


tan  ft 

x  sec2  at  +  tan2  ft 
sec  at 

x  sec2  a(  +  tan2  ft 
0 


tan  at 

x/  sec2  at  +  tan2  ft 

sin  a(  tan  ft 
x  sec2  at  +  tan2  Pi 

cos  at 


B4  Inverse  Transformations. 

There  is  no  need  to  list  the  inverse  transformations,  e.g.  from  wind  axes  to  body  axes 
M]nw.  since 


and  making  use  of  the  property  of  direction  cosines  that 


[A]1  =  [Ali- 

then 

[A]bw  =  [A]wbt 

i.e.  the  inverse  of  a  transformation  is  obtained  by  simply  taking  the  transpose  of  the  direct  matrix. 


B5. — Transformations  Between  Other  Systems. 

Transformations  between  the  stability-  and  wind-axes  systems  may  be  derived  from  the 
matrices  listed  above  as  follows.  We  wish  to  find  the  matrix  [A)sw  such  that 

As  =  [A]Sw  Aw 

Now 

Aw  —  [AIwhAb 

Hence 

A.s  -  [AJsw  [A]wb  Ab 
But  we  also  know  that 
As  -  [A]sb  An 

Hence 

(A|sw  [A|wb  —  [A]sb 
i.e. 

[A)sw  —  [AIsb  [A]wb  1 
[A]sr  [A)wiit 

Hence  in  terms  of  ott  and  (the  most  useful  pair  for  this  transformation) 


[A)sw 


cos  &  sin  /J  0 

sin  f)  cos  fl  0 

0  0  1 


Similar  transformations  for  other  *,  /S  pairs  may  be  derived  from  the  matrices  listed  above. 
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